
     Program PWSCF v.6.3MaX starts on 13Oct2018 at 10: 4: 3 

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
         "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);
          URL http://www.quantum-espresso.org", 
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/quote

     Serial multi-threaded version, running on    6 processor cores
     Waiting for input...
     Reading input from standard input

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
               file Bi.pz-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s)  5D renormalized
     Found symmetry operation: I + ( -0.3333  0.3333  0.3333)
     This is a supercell, fractional translations are disabled

     G-vector sticks info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Sum        2515    1261    361               252985    89447   13545



     bravais-lattice index     =            0
     lattice parameter (alat)  =      10.6547  a.u.
     unit-cell volume          =    2614.6950 (a.u.)^3
     number of atoms/cell      =           30
     number of atomic types    =            4
     number of electrons       =       294.00
     number of Kohn-Sham states=          176
     kinetic-energy cutoff     =      40.0000  Ry
     charge density cutoff     =     320.0000  Ry
     convergence threshold     =      1.0E-06
     mixing beta               =       0.4000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      =  SLA  PZ   NOGX NOGC ( 1  1  0  0 0 0)

     celldm(1)=  10.654661  celldm(2)=   0.000000  celldm(3)=   0.000000
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   1.000000   0.000000  -0.000000 )  
               a(2) = (  -0.500000   0.866025  -0.000000 )  
               a(3) = (   0.000000   0.000000   2.496155 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  1.000000  0.577350 -0.000000 )  
               b(2) = ( -0.000000  1.154701  0.000000 )  
               b(3) = (  0.000000  0.000000  0.400616 )  


     PseudoPot. # 1 for Bi read from file:
     /home/ngt/.burai/.pseudopot/Bi.pz-dn-rrkjus_psl.1.0.0.UPF
     MD5 check sum: 30ed9582656127c1bc49cf5897a59f42
     Pseudo is Ultrasoft + core correction, Zval = 15.0
     Generated using "atomic" code by A. Dal Corso v.6.2.2
     Using radial grid of 1283 points,  6 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
                l(5) =   2
                l(6) =   2
     Q(r) pseudized with 0 coefficients 


     PseudoPot. # 2 for O  read from file:
     /home/ngt/.burai/.pseudopot/o_lda_v1.2.uspp.F.UPF
     MD5 check sum: 194b818514d2a7dbb603bb78519a4f40
     Pseudo is Ultrasoft + core correction, Zval =  6.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  737 points,  5 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
                l(5) =   2
     Q(r) pseudized with  8 coefficients,  rinner =    0.900   0.900   0.900
                                                       0.900   0.900

     PseudoPot. # 3 for Fe read from file:
     /home/ngt/.burai/.pseudopot/fe_lda_v1.5.uspp.F.UPF
     MD5 check sum: 70aee299febfc78ac0df086b2000f8f7
     Pseudo is Ultrasoft + core correction, Zval = 16.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  861 points,  6 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
                l(5) =   2
                l(6) =   2
     Q(r) pseudized with  8 coefficients,  rinner =    1.000   1.000   1.000
                                                       1.000   1.000

     PseudoPot. # 4 for Fe read from file:
     /home/ngt/.burai/.pseudopot/fe_lda_v1.5.uspp.F.UPF
     MD5 check sum: 70aee299febfc78ac0df086b2000f8f7
     Pseudo is Ultrasoft + core correction, Zval = 16.0
     Generated by new atomic code, or converted to UPF format
     Using radial grid of  861 points,  6 beta functions with: 
                l(1) =   0
                l(2) =   0
                l(3) =   1
                l(4) =   1
                l(5) =   2
                l(6) =   2
     Q(r) pseudized with  8 coefficients,  rinner =    1.000   1.000   1.000
                                                       1.000   1.000

     atomic species   valence    mass     pseudopotential
        Bi            15.00   208.98038     Bi( 1.00)
        O              6.00    15.99940     O ( 1.00)
        Fe1           16.00    55.84500     Fe( 1.00)
        Fe2           16.00    55.84500     Fe( 1.00)

     Starting magnetic structure 
     atomic species   magnetization
        Bi           0.000
        O            0.000
        Fe1          0.800
        Fe2         -0.800


     Simplified LDA+U calculation (l_max = 2) with parameters (eV):
     atomic species    L          U    alpha       J0     beta
        Fe1            2     2.4346   0.0000   0.0000   0.0000
        Fe2            2     2.4346   0.0000   0.0000   0.0000



      3 Sym. Ops. (no inversion) found


                                    s                        frac. trans.

      isym =  1     identity                                     

 cryst.   s( 1) = (     1          0          0      )
                  (     0          1          0      )
                  (     0          0          1      )

 cart.    s( 1) = (  1.0000000  0.0000000  0.0000000 )
                  (  0.0000000  1.0000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


      isym =  2     120 deg rotation - cryst. axis [0,0,1]       

 cryst.   s( 2) = (     0          1          0      )
                  (    -1         -1          0      )
                  (     0          0          1      )

 cart.    s( 2) = ( -0.5000000 -0.8660254  0.0000000 )
                  (  0.8660254 -0.5000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


      isym =  3     120 deg rotation - cryst. axis [0,0,-1]      

 cryst.   s( 3) = (    -1         -1          0      )
                  (     1          0          0      )
                  (     0          0          1      )

 cart.    s( 3) = ( -0.5000000  0.8660254  0.0000000 )
                  ( -0.8660254 -0.5000000  0.0000000 )
                  (  0.0000000  0.0000000  1.0000000 )


   Cartesian axes

     site n.     atom                  positions (alat units)
         1           Fe1 tau(   1) = (   0.0000000   0.0000000   0.6985776  )
         2           Fe2 tau(   2) = (   0.5000005   0.2886748   0.2825564  )
         3           Fe1 tau(   3) = (   0.5000005   0.2886748   1.5306360  )
         4           Fe2 tau(   4) = (  -0.0000005   0.5773505   1.1145958  )
         5           Fe1 tau(   5) = (  -0.0000005   0.5773505   2.3626695  )
         6           Fe2 tau(   6) = (   0.0000000   0.0000000   1.9466681  )
         7           Bi  tau(   7) = (  -0.0000005   0.5773505   0.4200321  )
         8           Bi  tau(   8) = (   0.0000000   0.0000000   0.0040159  )
         9           Bi  tau(   9) = (   0.0000000   0.0000000   1.2520984  )
        10           Bi  tau(  10) = (   0.5000005   0.2886748   0.8360602  )
        11           Bi  tau(  11) = (   0.5000005   0.2886748   2.0841436  )
        12           Bi  tau(  12) = (  -0.0000005   0.5773505   1.6681117  )
        13           O   tau(  13) = (   0.6984987   0.0853327   0.5387765  )
        14           O   tau(  14) = (   0.5768504   0.5622513   0.5387765  )
        15           O   tau(  15) = (   0.2246509   0.2184414   0.5387765  )
        16           O   tau(  16) = (   0.2753499   0.5071158   0.1227567  )
        17           O   tau(  17) = (  -0.1984997   0.3740075   0.1227567  )
        18           O   tau(  18) = (  -0.0768501   0.8509274   0.1227567  )
        19           O   tau(  19) = (   0.1984987   0.3740076   1.3708274  )
        20           O   tau(  20) = (   0.0768506   0.8509265   1.3708274  )
        21           O   tau(  21) = (  -0.2753492   0.5071167   1.3708274  )
        22           O   tau(  22) = (  -0.2246502   0.7957908   0.9547978  )
        23           O   tau(  23) = (   0.3015002   0.6626827   0.9547978  )
        24           O   tau(  24) = (  -0.0768500   0.2735772   0.9547978  )
        25           O   tau(  25) = (  -0.3015049   0.6626846   2.2028821  )
        26           O   tau(  26) = (   0.0768507   0.2735721   2.2028821  )
        27           O   tau(  27) = (   0.2246542   0.7957940   2.2028821  )
        28           O   tau(  28) = (   0.7753507   0.2184423   1.7868570  )
        29           O   tau(  29) = (   0.3015012   0.0853309   1.7868570  )
        30           O   tau(  30) = (   0.4231481   0.5622522   1.7868570  )

   Crystallographic axes

     site n.     atom                  positions (cryst. coord.)
         1           Fe1 tau(   1) = (  0.0000000  0.0000000  0.2798615  )
         2           Fe2 tau(   2) = (  0.6666670  0.3333330  0.1131967  )
         3           Fe1 tau(   3) = (  0.6666670  0.3333330  0.6131975  )
         4           Fe2 tau(   4) = (  0.3333330  0.6666670  0.4465251  )
         5           Fe1 tau(   5) = (  0.3333330  0.6666670  0.9465236  )
         6           Fe2 tau(   6) = (  0.0000000  0.0000000  0.7798667  )
         7           Bi  tau(   7) = (  0.3333330  0.6666670  0.1682716  )
         8           Bi  tau(   8) = (  0.0000000  0.0000000  0.0016088  )
         9           Bi  tau(   9) = (  0.0000000  0.0000000  0.5016109  )
        10           Bi  tau(  10) = (  0.6666670  0.3333330  0.3349393  )
        11           Bi  tau(  11) = (  0.6666670  0.3333330  0.8349416  )
        12           Bi  tau(  12) = (  0.3333330  0.6666670  0.6682725  )
        13           O   tau(  13) = (  0.7477656  0.0985337  0.2158426  )
        14           O   tau(  14) = (  0.9014663  0.6492319  0.2158426  )
        15           O   tau(  15) = (  0.3507681  0.2522344  0.2158426  )
        16           O   tau(  16) = (  0.5681333  0.5855669  0.0491783  )
        17           O   tau(  17) = (  0.0174336  0.4318667  0.0491783  )
        18           O   tau(  18) = (  0.4144331  0.9825664  0.0491783  )
        19           O   tau(  19) = (  0.4144321  0.4318668  0.5491756  )
        20           O   tau(  20) = (  0.5681332  0.9825653  0.5491756  )
        21           O   tau(  21) = (  0.0174347  0.5855679  0.5491756  )
        22           O   tau(  22) = (  0.2347999  0.9189001  0.3825074  )
        23           O   tau(  23) = (  0.6841002  0.7652001  0.3825074  )
        24           O   tau(  24) = (  0.0810999  0.3158998  0.3825074  )
        25           O   tau(  25) = (  0.0810962  0.7652023  0.8825102  )
        26           O   tau(  26) = (  0.2347977  0.3158939  0.8825102  )
        27           O   tau(  27) = (  0.6841061  0.9189038  0.8825102  )
        28           O   tau(  28) = (  0.9014684  0.2522354  0.7158438  )
        29           O   tau(  29) = (  0.3507670  0.0985316  0.7158438  )
        30           O   tau(  30) = (  0.7477646  0.6492330  0.7158438  )

     number of k points=     8  gaussian smearing, width (Ry)=  0.0100
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0370370
        k(    2) = (   0.0000000   0.0000000   0.1335387), wk =   0.0740741
        k(    3) = (   0.0000000   0.3849002   0.0000000), wk =   0.2222222
        k(    4) = (   0.0000000   0.3849002   0.1335387), wk =   0.2222222
        k(    5) = (   0.3333333   0.5773503   0.0000000), wk =   0.0740741
        k(    6) = (   0.3333333   0.5773503   0.1335387), wk =   0.0740741
        k(    7) = (  -0.3333333   0.1924501   0.1335387), wk =   0.2222222
        k(    8) = (   0.3333333  -0.5773503  -0.1335387), wk =   0.0740741

                       cryst. coord.
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0370370
        k(    2) = (   0.0000000   0.0000000   0.3333333), wk =   0.0740741
        k(    3) = (   0.0000000   0.3333333   0.0000000), wk =   0.2222222
        k(    4) = (   0.0000000   0.3333333   0.3333333), wk =   0.2222222
        k(    5) = (   0.3333333   0.3333333   0.0000000), wk =   0.0740741
        k(    6) = (   0.3333333   0.3333333   0.3333333), wk =   0.0740741
        k(    7) = (  -0.3333333   0.3333333   0.3333333), wk =   0.2222222
        k(    8) = (   0.3333333  -0.6666667  -0.3333333), wk =   0.0740741

     Dense  grid:   252985 G-vectors     FFT dimensions: (  64,  64, 160)

     Smooth grid:    89447 G-vectors     FFT dimensions: (  45,  45, 108)

     Dynamical RAM for                 wfc:      30.00 MB

     Dynamical RAM for     wfc (w. buffer):     509.96 MB

     Dynamical RAM for             U proj.:       5.11 MB

     Dynamical RAM for  U proj. (w. buff.):      86.92 MB

     Dynamical RAM for           str. fact:      15.44 MB

     Dynamical RAM for           local pot:       0.00 MB

     Dynamical RAM for          nlocal pot:      76.70 MB

     Dynamical RAM for                qrad:       9.73 MB

     Dynamical RAM for          rho,v,vnew:      53.16 MB

     Dynamical RAM for               rhoin:      17.72 MB

     Dynamical RAM for            rho*nmix:     123.53 MB

     Dynamical RAM for           G-vectors:      15.16 MB

     Dynamical RAM for          h,s,v(r/c):       1.42 MB

     Dynamical RAM for          <psi|beta>:       1.21 MB

     Dynamical RAM for      wfcinit/wfcrot:      70.49 MB

     Dynamical RAM for           addusdens:     791.35 MB

     Estimated static dynamical RAM per process >     957.62 MB

     Estimated max dynamical RAM per process >       1.71 GB
     Generating pointlists ...
     new r_m :   0.1696 (alat units)  1.8074 (a.u.) for type    1
     new r_m :   0.1451 (alat units)  1.5459 (a.u.) for type    2
     new r_m :   0.1451 (alat units)  1.5459 (a.u.) for type    3
     new r_m :   0.1451 (alat units)  1.5459 (a.u.) for type    4

     Check: negative/imaginary core charge=   -0.000012    0.000000

     Initial potential from superposition of free atoms

     starting charge  287.99607, renormalised to  294.00000
 --- in v_hubbard ---
Hubbard energy    0.4295
 -------
     Number of +U iterations with fixed ns =  0
     Starting occupations:
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   5.00000  1.00000  6.00000
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
   spin  2
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
atomic mag. moment =     4.000000
atom    2   Tr[ns(na)] (up, down, total) =   1.00000  5.00000  6.00000
   spin  1
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
atomic mag. moment =    -4.000000
atom    3   Tr[ns(na)] (up, down, total) =   5.00000  1.00000  6.00000
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
   spin  2
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
atomic mag. moment =     4.000000
atom    4   Tr[ns(na)] (up, down, total) =   1.00000  5.00000  6.00000
   spin  1
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
atomic mag. moment =    -4.000000
atom    5   Tr[ns(na)] (up, down, total) =   5.00000  1.00000  6.00000
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
   spin  2
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
atomic mag. moment =     4.000000
atom    6   Tr[ns(na)] (up, down, total) =   1.00000  5.00000  6.00000
   spin  1
    eigenvalues: 
  0.200  0.200  0.200  0.200  0.200
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.200  0.000  0.000  0.000
  0.000  0.000  0.200  0.000  0.000
  0.000  0.000  0.000  0.200  0.000
  0.000  0.000  0.000  0.000  0.200
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.000  1.000
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.000  0.000
  0.000  0.000  1.000  0.000  0.000
  0.000  0.000  0.000  1.000  0.000
  0.000  0.000  0.000  0.000  1.000
atomic mag. moment =    -4.000000
N of occupied +U levels =   36.000000
 --- exit write_ns ---
 Atomic wfc used for LDA+U Projector are NOT orthogonalized
     Starting wfcs are  204 randomized atomic wfcs

     total cpu time spent up to now is       25.1 secs

     Self-consistent Calculation

     iteration #  1     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  1.00E-02,  avg # of iterations =  2.3
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   5.00230  0.29549  5.29779
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.039  0.258  0.000  0.349  0.354
  0.258  0.039  0.000  0.354  0.349
  0.572  0.131  0.000  0.297  0.000
  0.131  0.572  0.000  0.000  0.297
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000  0.000 -0.000
  0.000  0.000  1.001 -0.000 -0.000
  0.000  0.000 -0.000  1.000 -0.000
  0.000 -0.000 -0.000 -0.000  1.000
   spin  2
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.315  0.000  0.146  0.453
  0.315  0.087  0.000  0.453  0.146
  0.559  0.040  0.000  0.381  0.020
  0.040  0.559  0.000  0.020  0.381
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000 -0.022  0.024
  0.000  0.000  0.072  0.024  0.022
  0.000 -0.022  0.024  0.059 -0.000
  0.000  0.024  0.022 -0.000  0.059
atomic mag. moment =     4.706809
atom    2   Tr[ns(na)] (up, down, total) =   0.29544  5.00229  5.29773
   spin  1
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.356  0.044  0.000  0.594  0.005
  0.044  0.356  0.000  0.005  0.594
  0.468  0.131  0.000  0.221  0.179
  0.131  0.468  0.000  0.179  0.221
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000  0.022  0.024
  0.000  0.000  0.072  0.024 -0.022
  0.000  0.022  0.024  0.059 -0.000
  0.000  0.024 -0.022 -0.000  0.059
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.287  0.010  0.000  0.700  0.003
  0.010  0.287  0.000  0.003  0.700
  0.493  0.209  0.000  0.174  0.124
  0.209  0.493  0.000  0.124  0.174
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000 -0.000 -0.000
  0.000  0.000  1.001 -0.000  0.000
  0.000 -0.000 -0.000  1.000 -0.000
  0.000 -0.000  0.000 -0.000  1.000
atomic mag. moment =    -4.706850
atom    3   Tr[ns(na)] (up, down, total) =   5.00228  0.29520  5.29748
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.039  0.256  0.000  0.352  0.353
  0.256  0.039  0.000  0.353  0.352
  0.573  0.132  0.000  0.295  0.000
  0.132  0.573  0.000  0.000  0.295
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000  0.000 -0.000
  0.000  0.000  1.001 -0.000 -0.000
  0.000  0.000 -0.000  1.000 -0.000
  0.000 -0.000 -0.000 -0.000  1.000
   spin  2
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.086  0.315  0.000  0.145  0.453
  0.315  0.086  0.000  0.453  0.145
  0.559  0.040  0.000  0.381  0.020
  0.040  0.559  0.000  0.020  0.381
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000 -0.022  0.024
  0.000  0.000  0.072  0.024  0.022
  0.000 -0.022  0.024  0.059 -0.000
  0.000  0.024  0.022 -0.000  0.059
atomic mag. moment =     4.707081
atom    4   Tr[ns(na)] (up, down, total) =   0.29543  5.00228  5.29771
   spin  1
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.358  0.044  0.000  0.593  0.005
  0.044  0.358  0.000  0.005  0.593
  0.466  0.131  0.000  0.222  0.180
  0.131  0.466  0.000  0.180  0.222
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000  0.022  0.024
  0.000  0.000  0.072  0.024 -0.022
  0.000  0.022  0.024  0.059 -0.000
  0.000  0.024 -0.022 -0.000  0.059
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.287  0.010  0.000  0.700  0.003
  0.010  0.287  0.000  0.003  0.700
  0.494  0.209  0.000  0.174  0.123
  0.209  0.494  0.000  0.124  0.174
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000 -0.000 -0.000
  0.000  0.000  1.001 -0.000  0.000
  0.000 -0.000 -0.000  1.000 -0.000
  0.000 -0.000  0.000 -0.000  1.000
atomic mag. moment =    -4.706853
atom    5   Tr[ns(na)] (up, down, total) =   5.00228  0.29524  5.29753
   spin  1
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.038  0.258  0.000  0.349  0.354
  0.258  0.038  0.000  0.355  0.349
  0.572  0.131  0.000  0.296  0.000
  0.131  0.572  0.000  0.000  0.296
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000  0.000 -0.000
  0.000  0.000  1.001 -0.000 -0.000
  0.000  0.000 -0.000  1.000 -0.000
  0.000 -0.000 -0.000 -0.000  1.000
   spin  2
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.314  0.000  0.146  0.452
  0.314  0.087  0.000  0.452  0.146
  0.559  0.039  0.000  0.381  0.020
  0.039  0.559  0.000  0.020  0.381
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000 -0.022  0.024
  0.000  0.000  0.072  0.024  0.022
  0.000 -0.022  0.024  0.059 -0.000
  0.000  0.024  0.022 -0.000  0.059
atomic mag. moment =     4.707039
atom    6   Tr[ns(na)] (up, down, total) =   0.29507  5.00229  5.29736
   spin  1
    eigenvalues: 
  0.032  0.032  0.033  0.099  0.099
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.356  0.044  0.000  0.594  0.005
  0.044  0.356  0.000  0.005  0.594
  0.468  0.131  0.000  0.221  0.179
  0.131  0.468  0.000  0.179  0.221
    occupations:
  0.033  0.000  0.000  0.000  0.000
  0.000  0.072  0.000  0.022  0.024
  0.000  0.000  0.072  0.024 -0.022
  0.000  0.022  0.024  0.059 -0.000
  0.000  0.024 -0.022 -0.000  0.059
   spin  2
    eigenvalues: 
  1.000  1.000  1.000  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.286  0.010  0.000  0.701  0.003
  0.010  0.286  0.000  0.003  0.701
  0.495  0.209  0.000  0.173  0.123
  0.209  0.494  0.000  0.123  0.173
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  1.001  0.000 -0.000 -0.000
  0.000  0.000  1.001 -0.000  0.000
  0.000 -0.000 -0.000  1.000 -0.000
  0.000 -0.000  0.000 -0.000  1.000
atomic mag. moment =    -4.707225
N of occupied +U levels =   31.785603
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.1451
 -------

     total cpu time spent up to now is      205.0 secs

     total energy              =   -3007.05911504 Ry
     Harris-Foulkes estimate   =   -3007.76635393 Ry
     estimated scf accuracy    <      18.97995378 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    30.12 Bohr mag/cell

     iteration #  2     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  6.46E-03,  avg # of iterations =  3.4
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   5.00030  1.04097  6.04127
   spin  1
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.155  0.431  0.000  0.107  0.308
  0.431  0.155  0.000  0.308  0.107
  0.391  0.024  0.000  0.550  0.036
  0.024  0.391  0.000  0.036  0.550
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.001 -0.001
  0.000  0.000  1.000 -0.001 -0.001
  0.000  0.001 -0.001  1.000 -0.000
  0.000 -0.001 -0.001 -0.000  1.000
   spin  2
    eigenvalues: 
  0.185  0.199  0.199  0.230  0.230
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.043  0.396  0.318  0.243
  0.000  0.396  0.043  0.243  0.318
  0.000  0.501  0.060  0.183  0.256
  0.000  0.060  0.501  0.256  0.183
    occupations:
  0.185  0.000  0.000  0.000  0.000
  0.000  0.216  0.000 -0.000  0.015
  0.000  0.000  0.216  0.015  0.000
  0.000 -0.000  0.015  0.212 -0.000
  0.000  0.015  0.000 -0.000  0.212
atomic mag. moment =     3.959333
atom    2   Tr[ns(na)] (up, down, total) =   1.03647  5.00030  6.03677
   spin  1
    eigenvalues: 
  0.186  0.199  0.199  0.226  0.226
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.066  0.357  0.366  0.211
  0.000  0.357  0.066  0.211  0.366
  0.000  0.499  0.078  0.166  0.257
  0.000  0.078  0.499  0.257  0.166
    occupations:
  0.186  0.000  0.000  0.000  0.000
  0.000  0.215  0.000  0.000  0.013
  0.000  0.000  0.215  0.013 -0.000
  0.000  0.000  0.013  0.210 -0.000
  0.000  0.013 -0.000 -0.000  0.210
   spin  2
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.581  0.004  0.000  0.378  0.037
  0.004  0.581  0.000  0.037  0.378
  0.239  0.176  0.000  0.458  0.127
  0.176  0.239  0.000  0.127  0.458
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000 -0.001 -0.001
  0.000  0.000  1.000 -0.001  0.001
  0.000 -0.001 -0.001  1.000 -0.000
  0.000 -0.001  0.001 -0.000  1.000
atomic mag. moment =    -3.963838
atom    3   Tr[ns(na)] (up, down, total) =   5.00031  1.04112  6.04142
   spin  1
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.154  0.431  0.000  0.107  0.308
  0.431  0.154  0.000  0.308  0.107
  0.391  0.024  0.000  0.549  0.036
  0.024  0.391  0.000  0.036  0.550
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.001 -0.001
  0.000  0.000  1.000 -0.001 -0.001
  0.000  0.001 -0.001  1.000 -0.000
  0.000 -0.001 -0.001 -0.000  1.000
   spin  2
    eigenvalues: 
  0.186  0.199  0.199  0.228  0.228
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.036  0.439  0.258  0.267
  0.000  0.439  0.036  0.267  0.258
  0.000  0.474  0.050  0.225  0.250
  0.000  0.050  0.474  0.250  0.225
    occupations:
  0.186  0.000  0.000  0.000  0.000
  0.000  0.214  0.000 -0.001  0.015
  0.000  0.000  0.214  0.015  0.001
  0.000 -0.001  0.015  0.213 -0.000
  0.000  0.015  0.001 -0.000  0.213
atomic mag. moment =     3.959190
atom    4   Tr[ns(na)] (up, down, total) =   1.04197  5.00029  6.04227
   spin  1
    eigenvalues: 
  0.186  0.198  0.198  0.230  0.230
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.054  0.397  0.311  0.238
  0.000  0.397  0.054  0.238  0.311
  0.000  0.483  0.066  0.196  0.255
  0.000  0.066  0.483  0.255  0.196
    occupations:
  0.186  0.000  0.000  0.000  0.000
  0.000  0.215  0.000 -0.000  0.016
  0.000  0.000  0.215  0.016  0.000
  0.000 -0.000  0.016  0.212 -0.000
  0.000  0.016  0.000 -0.000  0.212
   spin  2
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.583  0.004  0.000  0.375  0.038
  0.004  0.584  0.000  0.038  0.375
  0.237  0.176  0.000  0.460  0.127
  0.176  0.237  0.000  0.127  0.460
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000 -0.001 -0.001
  0.000  0.000  1.000 -0.001  0.001
  0.000 -0.001 -0.001  1.000 -0.000
  0.000 -0.001  0.001 -0.000  1.000
atomic mag. moment =    -3.958324
atom    5   Tr[ns(na)] (up, down, total) =   5.00031  1.03632  6.03663
   spin  1
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.155  0.431  0.000  0.106  0.308
  0.431  0.155  0.000  0.308  0.106
  0.390  0.024  0.000  0.550  0.036
  0.024  0.390  0.000  0.036  0.550
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000  0.001 -0.001
  0.000  0.000  1.000 -0.001 -0.001
  0.000  0.001 -0.001  1.000 -0.000
  0.000 -0.001 -0.001 -0.000  1.000
   spin  2
    eigenvalues: 
  0.185  0.199  0.199  0.227  0.227
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.034  0.410  0.292  0.263
  0.000  0.410  0.034  0.263  0.292
  0.000  0.499  0.056  0.192  0.252
  0.000  0.056  0.499  0.252  0.192
    occupations:
  0.185  0.000  0.000  0.000  0.000
  0.000  0.214  0.000 -0.001  0.014
  0.000  0.000  0.214  0.014  0.001
  0.000 -0.001  0.014  0.211 -0.000
  0.000  0.014  0.001 -0.000  0.211
atomic mag. moment =     3.963988
atom    6   Tr[ns(na)] (up, down, total) =   1.04033  5.00030  6.04063
   spin  1
    eigenvalues: 
  0.187  0.199  0.199  0.228  0.228
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.053  0.427  0.267  0.252
  0.000  0.427  0.053  0.252  0.267
  0.000  0.459  0.060  0.229  0.251
  0.000  0.060  0.459  0.251  0.229
    occupations:
  0.187  0.000  0.000  0.000  0.000
  0.000  0.214  0.000 -0.000  0.015
  0.000  0.000  0.214  0.015  0.000
  0.000 -0.000  0.015  0.213 -0.000
  0.000  0.015  0.000 -0.000  0.213
   spin  2
    eigenvalues: 
  0.999  0.999  1.001  1.001  1.001
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.582  0.004  0.000  0.377  0.037
  0.004  0.582  0.000  0.037  0.377
  0.238  0.176  0.000  0.459  0.127
  0.176  0.238  0.000  0.127  0.459
    occupations:
  1.001  0.000  0.000  0.000  0.000
  0.000  1.000  0.000 -0.001 -0.001
  0.000  0.000  1.000 -0.001  0.001
  0.000 -0.001 -0.001  1.000 -0.000
  0.000 -0.001  0.001 -0.000  1.000
atomic mag. moment =    -3.959968
N of occupied +U levels =   36.238983
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.2471
 -------

     total cpu time spent up to now is      398.1 secs

     total energy              =   -3010.79147549 Ry
     Harris-Foulkes estimate   =   -3010.71765734 Ry
     estimated scf accuracy    <       7.74098239 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    27.06 Bohr mag/cell

     iteration #  3     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  2.63E-03,  avg # of iterations =  3.3
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.99604  0.72449  5.72053
   spin  1
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.063  0.472  0.000  0.103  0.362
  0.472  0.063  0.000  0.362  0.103
  0.412  0.054  0.000  0.512  0.022
  0.054  0.412  0.000  0.022  0.512
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000  0.000 -0.001
  0.000  0.000  0.999 -0.001 -0.000
  0.000  0.000 -0.001  0.999 -0.000
  0.000 -0.001 -0.000 -0.000  0.999
   spin  2
    eigenvalues: 
  0.097  0.099  0.099  0.215  0.215
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.075  0.331  0.121  0.472
  0.000  0.331  0.075  0.472  0.121
  0.000  0.551  0.043  0.382  0.024
  0.000  0.043  0.551  0.024  0.382
    occupations:
  0.097  0.000  0.000  0.000  0.000
  0.000  0.168  0.000 -0.037  0.043
  0.000  0.000  0.168  0.043  0.037
  0.000 -0.037  0.043  0.146 -0.000
  0.000  0.043  0.037 -0.000  0.146
atomic mag. moment =     4.271542
atom    2   Tr[ns(na)] (up, down, total) =   0.72475  4.99610  5.72085
   spin  1
    eigenvalues: 
  0.098  0.099  0.099  0.214  0.214
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.359  0.046  0.588  0.007
  0.000  0.046  0.359  0.007  0.588
  0.000  0.456  0.139  0.221  0.184
  0.000  0.139  0.456  0.184  0.221
    occupations:
  0.098  0.000  0.000  0.000  0.000
  0.000  0.168  0.000  0.037  0.042
  0.000  0.000  0.168  0.042 -0.037
  0.000  0.037  0.042  0.146 -0.000
  0.000  0.042 -0.037 -0.000  0.146
   spin  2
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.521  0.014  0.000  0.405  0.060
  0.014  0.521  0.000  0.060  0.405
  0.268  0.197  0.000  0.410  0.125
  0.197  0.267  0.000  0.125  0.410
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000 -0.000 -0.001
  0.000  0.000  0.999 -0.001  0.000
  0.000 -0.000 -0.001  0.999 -0.000
  0.000 -0.001  0.000 -0.000  0.999
atomic mag. moment =    -4.271349
atom    3   Tr[ns(na)] (up, down, total) =   4.99602  0.72291  5.71893
   spin  1
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.061  0.472  0.000  0.104  0.363
  0.471  0.061  0.000  0.363  0.104
  0.412  0.055  0.000  0.511  0.022
  0.055  0.412  0.000  0.022  0.511
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000  0.000 -0.001
  0.000  0.000  0.999 -0.001 -0.000
  0.000  0.000 -0.001  0.999 -0.000
  0.000 -0.001 -0.000 -0.000  0.999
   spin  2
    eigenvalues: 
  0.097  0.099  0.099  0.214  0.214
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.076  0.333  0.119  0.472
  0.000  0.333  0.076  0.472  0.119
  0.000  0.549  0.042  0.383  0.025
  0.000  0.042  0.549  0.025  0.383
    occupations:
  0.097  0.000  0.000  0.000  0.000
  0.000  0.167  0.000 -0.037  0.043
  0.000  0.000  0.167  0.043  0.037
  0.000 -0.037  0.043  0.146 -0.000
  0.000  0.043  0.037 -0.000  0.146
atomic mag. moment =     4.273116
atom    4   Tr[ns(na)] (up, down, total) =   0.72506  4.99602  5.72108
   spin  1
    eigenvalues: 
  0.097  0.099  0.099  0.215  0.215
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.358  0.048  0.586  0.008
  0.000  0.048  0.358  0.008  0.586
  0.000  0.455  0.139  0.221  0.186
  0.000  0.139  0.455  0.186  0.221
    occupations:
  0.097  0.000  0.000  0.000  0.000
  0.000  0.168  0.000  0.038  0.043
  0.000  0.000  0.168  0.043 -0.038
  0.000  0.038  0.043  0.146 -0.000
  0.000  0.043 -0.038 -0.000  0.146
   spin  2
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.519  0.014  0.000  0.406  0.060
  0.014  0.519  0.000  0.060  0.406
  0.269  0.198  0.000  0.409  0.125
  0.198  0.269  0.000  0.125  0.409
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000 -0.000 -0.000
  0.000  0.000  0.999 -0.000  0.000
  0.000 -0.000 -0.000  0.999 -0.000
  0.000 -0.000  0.000 -0.000  0.999
atomic mag. moment =    -4.270963
atom    5   Tr[ns(na)] (up, down, total) =   4.99612  0.72597  5.72209
   spin  1
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.061  0.476  0.000  0.100  0.363
  0.475  0.061  0.000  0.363  0.100
  0.410  0.054  0.000  0.513  0.023
  0.054  0.410  0.000  0.023  0.514
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000  0.000 -0.001
  0.000  0.000  0.999 -0.001 -0.000
  0.000  0.000 -0.001  0.999 -0.000
  0.000 -0.001 -0.000 -0.000  0.999
   spin  2
    eigenvalues: 
  0.098  0.099  0.099  0.215  0.215
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.076  0.330  0.123  0.470
  0.000  0.330  0.076  0.470  0.123
  0.000  0.551  0.043  0.382  0.024
  0.000  0.043  0.551  0.024  0.382
    occupations:
  0.098  0.000  0.000  0.000  0.000
  0.000  0.168  0.000 -0.037  0.043
  0.000  0.000  0.168  0.043  0.037
  0.000 -0.037  0.043  0.146 -0.000
  0.000  0.043  0.037 -0.000  0.146
atomic mag. moment =     4.270151
atom    6   Tr[ns(na)] (up, down, total) =   0.72372  4.99603  5.71975
   spin  1
    eigenvalues: 
  0.097  0.099  0.099  0.214  0.214
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.360  0.048  0.584  0.008
  0.000  0.048  0.360  0.008  0.584
  0.000  0.454  0.138  0.222  0.186
  0.000  0.138  0.454  0.186  0.222
    occupations:
  0.097  0.000  0.000  0.000  0.000
  0.000  0.167  0.000  0.037  0.043
  0.000  0.000  0.167  0.043 -0.037
  0.000  0.037  0.043  0.146 -0.000
  0.000  0.043 -0.037 -0.000  0.146
   spin  2
    eigenvalues: 
  0.998  0.998  1.000  1.000  1.000
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.518  0.015  0.000  0.404  0.063
  0.015  0.518  0.000  0.063  0.404
  0.268  0.199  0.000  0.409  0.125
  0.199  0.268  0.000  0.125  0.409
    occupations:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.999  0.000 -0.000 -0.001
  0.000  0.000  0.999 -0.001  0.000
  0.000 -0.000 -0.001  0.999 -0.000
  0.000 -0.001  0.000 -0.000  0.999
atomic mag. moment =    -4.272316
N of occupied +U levels =   34.323228
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.3321
 -------

     total cpu time spent up to now is      584.8 secs

     total energy              =   -3012.75737796 Ry
     Harris-Foulkes estimate   =   -3011.65874316 Ry
     estimated scf accuracy    <       2.17014439 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    27.82 Bohr mag/cell

     iteration #  4     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  7.38E-04,  avg # of iterations =  3.5
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.98661  1.22291  6.20952
   spin  1
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.157  0.442  0.000  0.094  0.307
  0.442  0.157  0.000  0.307  0.094
  0.381  0.019  0.000  0.561  0.038
  0.019  0.381  0.000  0.038  0.561
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.997  0.000  0.001 -0.002
  0.000  0.000  0.997 -0.002 -0.001
  0.000  0.001 -0.002  0.997 -0.000
  0.000 -0.002 -0.001 -0.000  0.997
   spin  2
    eigenvalues: 
  0.198  0.203  0.203  0.309  0.309
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.052  0.343  0.091  0.514
  0.000  0.343  0.052  0.514  0.091
  0.000  0.556  0.048  0.369  0.026
  0.000  0.048  0.556  0.026  0.369
    occupations:
  0.198  0.000  0.000  0.000  0.000
  0.000  0.267  0.000 -0.032  0.041
  0.000  0.000  0.267  0.041  0.032
  0.000 -0.032  0.041  0.245 -0.000
  0.000  0.041  0.032 -0.000  0.245
atomic mag. moment =     3.763697
atom    2   Tr[ns(na)] (up, down, total) =   1.21290  4.98649  6.19939
   spin  1
    eigenvalues: 
  0.198  0.202  0.202  0.306  0.306
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.338  0.054  0.597  0.010
  0.000  0.054  0.338  0.010  0.597
  0.000  0.463  0.144  0.205  0.187
  0.000  0.144  0.463  0.187  0.205
    occupations:
  0.198  0.000  0.000  0.000  0.000
  0.000  0.265  0.000  0.032  0.039
  0.000  0.000  0.265  0.039 -0.032
  0.000  0.032  0.039  0.243 -0.000
  0.000  0.039 -0.032 -0.000  0.243
   spin  2
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.596  0.005  0.000  0.355  0.044
  0.005  0.596  0.000  0.044  0.355
  0.222  0.177  0.000  0.474  0.127
  0.177  0.222  0.000  0.127  0.474
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.996  0.000 -0.001 -0.002
  0.000  0.000  0.996 -0.002  0.001
  0.000 -0.001 -0.002  0.997 -0.000
  0.000 -0.002  0.001 -0.000  0.997
atomic mag. moment =    -3.773584
atom    3   Tr[ns(na)] (up, down, total) =   4.98675  1.22317  6.20992
   spin  1
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.156  0.441  0.000  0.095  0.307
  0.441  0.156  0.000  0.308  0.095
  0.383  0.020  0.000  0.560  0.038
  0.020  0.383  0.000  0.038  0.560
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.997  0.000  0.001 -0.002
  0.000  0.000  0.997 -0.002 -0.001
  0.000  0.001 -0.002  0.997 -0.000
  0.000 -0.002 -0.001 -0.000  0.997
   spin  2
    eigenvalues: 
  0.200  0.205  0.205  0.307  0.307
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.056  0.348  0.086  0.510
  0.000  0.348  0.056  0.510  0.086
  0.000  0.550  0.046  0.374  0.030
  0.000  0.046  0.550  0.030  0.374
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.266  0.000 -0.031  0.039
  0.000  0.000  0.266  0.039  0.031
  0.000 -0.031  0.039  0.246 -0.000
  0.000  0.039  0.031 -0.000  0.246
atomic mag. moment =     3.763587
atom    4   Tr[ns(na)] (up, down, total) =   1.22263  4.98651  6.20915
   spin  1
    eigenvalues: 
  0.199  0.203  0.203  0.309  0.309
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.335  0.062  0.591  0.013
  0.000  0.062  0.335  0.013  0.591
  0.000  0.460  0.144  0.204  0.192
  0.000  0.144  0.460  0.192  0.204
    occupations:
  0.199  0.000  0.000  0.000  0.000
  0.000  0.267  0.000  0.032  0.041
  0.000  0.000  0.267  0.041 -0.032
  0.000  0.032  0.041  0.245 -0.000
  0.000  0.041 -0.032 -0.000  0.245
   spin  2
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.594  0.005  0.000  0.356  0.045
  0.005  0.594  0.000  0.045  0.356
  0.223  0.178  0.000  0.472  0.127
  0.178  0.223  0.000  0.127  0.472
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.996  0.000 -0.001 -0.002
  0.000  0.000  0.996 -0.002  0.001
  0.000 -0.001 -0.002  0.997 -0.000
  0.000 -0.002  0.001 -0.000  0.997
atomic mag. moment =    -3.763879
atom    5   Tr[ns(na)] (up, down, total) =   4.98643  1.20796  6.19439
   spin  1
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.155  0.444  0.000  0.094  0.308
  0.444  0.155  0.000  0.308  0.094
  0.382  0.020  0.000  0.561  0.038
  0.020  0.382  0.000  0.038  0.561
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.996  0.000  0.001 -0.002
  0.000  0.000  0.996 -0.002 -0.001
  0.000  0.001 -0.002  0.997 -0.000
  0.000 -0.002 -0.001 -0.000  0.997
   spin  2
    eigenvalues: 
  0.196  0.201  0.201  0.305  0.305
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.056  0.338  0.105  0.501
  0.000  0.338  0.056  0.501  0.105
  0.000  0.558  0.049  0.371  0.023
  0.000  0.049  0.558  0.023  0.371
    occupations:
  0.196  0.000  0.000  0.000  0.000
  0.000  0.264  0.000 -0.032  0.040
  0.000  0.000  0.264  0.040  0.032
  0.000 -0.032  0.040  0.242 -0.000
  0.000  0.040  0.032 -0.000  0.242
atomic mag. moment =     3.778467
atom    6   Tr[ns(na)] (up, down, total) =   1.22193  4.98679  6.20872
   spin  1
    eigenvalues: 
  0.200  0.204  0.204  0.307  0.307
    eigenvectors:
  1.000  0.000  0.000  0.000  0.000
  0.000  0.341  0.060  0.587  0.012
  0.000  0.060  0.341  0.012  0.587
  0.000  0.456  0.143  0.209  0.192
  0.000  0.143  0.456  0.192  0.209
    occupations:
  0.200  0.000  0.000  0.000  0.000
  0.000  0.266  0.000  0.031  0.040
  0.000  0.000  0.266  0.040 -0.031
  0.000  0.031  0.040  0.245 -0.000
  0.000  0.040 -0.031 -0.000  0.245
   spin  2
    eigenvalues: 
  0.995  0.995  0.999  0.999  0.999
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.592  0.004  0.000  0.362  0.042
  0.004  0.592  0.000  0.042  0.362
  0.227  0.177  0.000  0.469  0.127
  0.177  0.227  0.000  0.127  0.469
    occupations:
  0.999  0.000  0.000  0.000  0.000
  0.000  0.997  0.000 -0.001 -0.001
  0.000  0.000  0.997 -0.001  0.001
  0.000 -0.001 -0.001  0.997 -0.000
  0.000 -0.001  0.001 -0.000  0.997
atomic mag. moment =    -3.764861
N of occupied +U levels =   37.231095
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.3923
 -------

     total cpu time spent up to now is      790.5 secs

     total energy              =   -3012.99009215 Ry
     Harris-Foulkes estimate   =   -3013.03954507 Ry
     estimated scf accuracy    <       0.55558255 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    25.67 Bohr mag/cell

     iteration #  5     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     ethr =  1.89E-04,  avg # of iterations =  3.6
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96930  1.10750  6.07680
   spin  1
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.172  0.400  0.000  0.130  0.298
  0.400  0.172  0.000  0.298  0.130
  0.401  0.027  0.000  0.537  0.035
  0.027  0.401  0.000  0.035  0.537
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000  0.002 -0.002
  0.000  0.000  0.993 -0.002 -0.002
  0.000  0.002 -0.002  0.994 -0.000
  0.000 -0.002 -0.002 -0.000  0.994
   spin  2
    eigenvalues: 
  0.157  0.157  0.158  0.318  0.318
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.079  0.315  0.000  0.161  0.444
  0.315  0.079  0.000  0.444  0.161
  0.557  0.049  0.000  0.377  0.017
  0.049  0.557  0.000  0.017  0.377
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.254  0.000 -0.054  0.057
  0.000  0.000  0.254  0.057  0.054
  0.000 -0.054  0.057  0.220 -0.000
  0.000  0.057  0.054 -0.000  0.220
atomic mag. moment =     3.861805
atom    2   Tr[ns(na)] (up, down, total) =   1.11918  4.97022  6.08939
   spin  1
    eigenvalues: 
  0.159  0.159  0.161  0.320  0.320
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.358  0.036  0.000  0.602  0.004
  0.036  0.358  0.000  0.004  0.602
  0.464  0.141  0.000  0.218  0.176
  0.141  0.464  0.000  0.176  0.218
    occupations:
  0.161  0.000  0.000  0.000  0.000
  0.000  0.257  0.000  0.054  0.057
  0.000  0.000  0.257  0.057 -0.054
  0.000  0.054  0.057  0.223 -0.000
  0.000  0.057 -0.054 -0.000  0.223
   spin  2
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.571  0.002  0.000  0.400  0.028
  0.002  0.571  0.000  0.028  0.400
  0.255  0.172  0.000  0.446  0.126
  0.172  0.255  0.000  0.126  0.446
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000 -0.002 -0.002
  0.000  0.000  0.993 -0.002  0.002
  0.000 -0.002 -0.002  0.994 -0.000
  0.000 -0.002  0.002 -0.000  0.994
atomic mag. moment =    -3.851040
atom    3   Tr[ns(na)] (up, down, total) =   4.96936  1.10808  6.07744
   spin  1
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.170  0.401  0.000  0.130  0.299
  0.401  0.170  0.000  0.299  0.130
  0.401  0.027  0.000  0.537  0.035
  0.027  0.401  0.000  0.035  0.537
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000  0.002 -0.002
  0.000  0.000  0.993 -0.002 -0.002
  0.000  0.002 -0.002  0.994 -0.000
  0.000 -0.002 -0.002 -0.000  0.994
   spin  2
    eigenvalues: 
  0.157  0.157  0.159  0.317  0.317
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.080  0.317  0.000  0.158  0.445
  0.317  0.080  0.000  0.445  0.158
  0.555  0.049  0.000  0.379  0.018
  0.049  0.555  0.000  0.018  0.379
    occupations:
  0.159  0.000  0.000  0.000  0.000
  0.000  0.254  0.000 -0.053  0.057
  0.000  0.000  0.254  0.057  0.053
  0.000 -0.053  0.057  0.221 -0.000
  0.000  0.057  0.053 -0.000  0.221
atomic mag. moment =     3.861286
atom    4   Tr[ns(na)] (up, down, total) =   1.10990  4.96966  6.07957
   spin  1
    eigenvalues: 
  0.157  0.157  0.158  0.319  0.319
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.358  0.037  0.000  0.601  0.004
  0.037  0.358  0.000  0.004  0.601
  0.464  0.141  0.000  0.218  0.177
  0.141  0.464  0.000  0.177  0.218
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.255  0.000  0.054  0.058
  0.000  0.000  0.255  0.058 -0.054
  0.000  0.054  0.058  0.221 -0.000
  0.000  0.058 -0.054 -0.000  0.221
   spin  2
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.571  0.002  0.000  0.397  0.030
  0.002  0.571  0.000  0.030  0.397
  0.253  0.173  0.000  0.447  0.126
  0.173  0.253  0.000  0.126  0.447
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000 -0.002 -0.002
  0.000  0.000  0.993 -0.002  0.002
  0.000 -0.002 -0.002  0.994 -0.000
  0.000 -0.002  0.002 -0.000  0.994
atomic mag. moment =    -3.859763
atom    5   Tr[ns(na)] (up, down, total) =   4.97011  1.12341  6.09352
   spin  1
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.175  0.399  0.000  0.129  0.296
  0.399  0.175  0.000  0.296  0.129
  0.400  0.026  0.000  0.539  0.035
  0.026  0.400  0.000  0.035  0.539
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000  0.002 -0.002
  0.000  0.000  0.993 -0.002 -0.002
  0.000  0.002 -0.002  0.994 -0.000
  0.000 -0.002 -0.002 -0.000  0.994
   spin  2
    eigenvalues: 
  0.160  0.160  0.162  0.320  0.320
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.080  0.317  0.000  0.159  0.444
  0.317  0.080  0.000  0.444  0.159
  0.554  0.049  0.000  0.379  0.018
  0.049  0.554  0.000  0.018  0.379
    occupations:
  0.162  0.000  0.000  0.000  0.000
  0.000  0.257  0.000 -0.054  0.057
  0.000  0.000  0.257  0.057  0.054
  0.000 -0.054  0.057  0.224 -0.000
  0.000  0.057  0.054 -0.000  0.224
atomic mag. moment =     3.846705
atom    6   Tr[ns(na)] (up, down, total) =   1.11133  4.96948  6.08080
   spin  1
    eigenvalues: 
  0.158  0.158  0.159  0.318  0.318
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.360  0.037  0.000  0.598  0.004
  0.037  0.360  0.000  0.004  0.598
  0.461  0.141  0.000  0.219  0.178
  0.141  0.461  0.000  0.178  0.219
    occupations:
  0.159  0.000  0.000  0.000  0.000
  0.000  0.255  0.000  0.054  0.057
  0.000  0.000  0.255  0.057 -0.054
  0.000  0.054  0.057  0.222 -0.000
  0.000  0.057 -0.054 -0.000  0.222
   spin  2
    eigenvalues: 
  0.990  0.990  0.996  0.997  0.997
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.576  0.002  0.000  0.392  0.030
  0.002  0.576  0.000  0.030  0.392
  0.249  0.172  0.000  0.452  0.127
  0.172  0.249  0.000  0.127  0.452
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.993  0.000 -0.002 -0.002
  0.000  0.000  0.993 -0.002  0.002
  0.000 -0.002 -0.002  0.994 -0.000
  0.000 -0.002  0.002 -0.000  0.994
atomic mag. moment =    -3.858153
N of occupied +U levels =   36.497519
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4504
 -------

     total cpu time spent up to now is     1002.6 secs

     total energy              =   -3013.08013028 Ry
     Harris-Foulkes estimate   =   -3013.07904986 Ry
     estimated scf accuracy    <       0.04233333 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.79 Bohr mag/cell

     iteration #  6     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  4 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  3 eigenvalues not converged
     c_bands:  3 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  3 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     ethr =  1.44E-05,  avg # of iterations =  4.1
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96448  1.18089  6.14537
   spin  1
    eigenvalues: 
  0.988  0.988  0.996  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.188  0.399  0.000  0.122  0.291
  0.399  0.188  0.000  0.291  0.122
  0.391  0.022  0.000  0.549  0.038
  0.022  0.391  0.000  0.038  0.549
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.992  0.000  0.003 -0.003
  0.000  0.000  0.992 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.166  0.166  0.169  0.340  0.340
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.089  0.314  0.000  0.166  0.430
  0.314  0.089  0.000  0.430  0.166
  0.549  0.048  0.000  0.384  0.019
  0.048  0.549  0.000  0.019  0.384
    occupations:
  0.169  0.000  0.000  0.000  0.000
  0.000  0.270  0.000 -0.060  0.061
  0.000  0.000  0.270  0.061  0.060
  0.000 -0.060  0.061  0.236 -0.000
  0.000  0.061  0.060 -0.000  0.236
atomic mag. moment =     3.783586
atom    2   Tr[ns(na)] (up, down, total) =   1.12425  4.96045  6.08470
   spin  1
    eigenvalues: 
  0.155  0.155  0.158  0.328  0.328
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.367  0.032  0.000  0.597  0.003
  0.032  0.367  0.000  0.003  0.597
  0.459  0.141  0.000  0.225  0.175
  0.141  0.459  0.000  0.175  0.225
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.259  0.000  0.059  0.061
  0.000  0.000  0.259  0.061 -0.059
  0.000  0.059  0.061  0.224 -0.000
  0.000  0.061 -0.059 -0.000  0.224
   spin  2
    eigenvalues: 
  0.987  0.987  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.588  0.002  0.000  0.381  0.029
  0.002  0.588  0.000  0.029  0.381
  0.241  0.170  0.000  0.463  0.127
  0.170  0.241  0.000  0.127  0.463
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.836197
atom    3   Tr[ns(na)] (up, down, total) =   4.96432  1.17679  6.14111
   spin  1
    eigenvalues: 
  0.988  0.988  0.996  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.183  0.394  0.000  0.131  0.293
  0.394  0.183  0.000  0.293  0.131
  0.399  0.025  0.000  0.540  0.036
  0.025  0.399  0.000  0.036  0.540
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.992  0.000  0.003 -0.003
  0.000  0.000  0.992 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.165  0.165  0.169  0.339  0.339
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.309  0.000  0.177  0.427
  0.309  0.087  0.000  0.427  0.177
  0.554  0.050  0.000  0.379  0.017
  0.050  0.554  0.000  0.017  0.379
    occupations:
  0.169  0.000  0.000  0.000  0.000
  0.000  0.270  0.000 -0.060  0.061
  0.000  0.000  0.270  0.061  0.060
  0.000 -0.060  0.061  0.234 -0.000
  0.000  0.061  0.060 -0.000  0.234
atomic mag. moment =     3.787527
atom    4   Tr[ns(na)] (up, down, total) =   1.16785  4.96305  6.13090
   spin  1
    eigenvalues: 
  0.163  0.163  0.166  0.337  0.337
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.367  0.032  0.000  0.597  0.003
  0.032  0.367  0.000  0.003  0.597
  0.459  0.142  0.000  0.225  0.175
  0.142  0.459  0.000  0.175  0.225
    occupations:
  0.166  0.000  0.000  0.000  0.000
  0.000  0.268  0.000  0.060  0.061
  0.000  0.000  0.268  0.061 -0.060
  0.000  0.060  0.061  0.233 -0.000
  0.000  0.061 -0.060 -0.000  0.233
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.580  0.001  0.000  0.392  0.026
  0.001  0.580  0.000  0.026  0.392
  0.249  0.169  0.000  0.455  0.127
  0.169  0.249  0.000  0.127  0.455
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.795204
atom    5   Tr[ns(na)] (up, down, total) =   4.95973  1.10595  6.06568
   spin  1
    eigenvalues: 
  0.987  0.987  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.176  0.395  0.000  0.134  0.296
  0.395  0.176  0.000  0.296  0.134
  0.402  0.027  0.000  0.536  0.035
  0.027  0.402  0.000  0.035  0.536
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.151  0.151  0.154  0.325  0.325
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.084  0.307  0.000  0.181  0.428
  0.307  0.084  0.000  0.428  0.181
  0.558  0.051  0.000  0.376  0.015
  0.051  0.558  0.000  0.015  0.376
    occupations:
  0.154  0.000  0.000  0.000  0.000
  0.000  0.257  0.000 -0.059  0.061
  0.000  0.000  0.257  0.061  0.059
  0.000 -0.059  0.061  0.219 -0.000
  0.000  0.061  0.059 -0.000  0.219
atomic mag. moment =     3.853782
atom    6   Tr[ns(na)] (up, down, total) =   1.16545  4.96452  6.12997
   spin  1
    eigenvalues: 
  0.162  0.162  0.166  0.338  0.338
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.358  0.031  0.000  0.608  0.002
  0.031  0.358  0.000  0.002  0.608
  0.469  0.142  0.000  0.218  0.172
  0.142  0.469  0.000  0.172  0.218
    occupations:
  0.166  0.000  0.000  0.000  0.000
  0.000  0.269  0.000  0.060  0.061
  0.000  0.000  0.269  0.061 -0.060
  0.000  0.060  0.061  0.230 -0.000
  0.000  0.061 -0.060 -0.000  0.230
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.562  0.001  0.000  0.415  0.022
  0.001  0.562  0.000  0.022  0.415
  0.267  0.170  0.000  0.437  0.126
  0.170  0.267  0.000  0.126  0.437
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000 -0.003 -0.003
  0.000  0.000  0.992 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.799070
N of occupied +U levels =   36.697727
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4589
 -------

     total cpu time spent up to now is     1260.4 secs

     total energy              =   -3013.08792298 Ry
     Harris-Foulkes estimate   =   -3013.09133005 Ry
     estimated scf accuracy    <       0.01535708 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.58 Bohr mag/cell

     iteration #  7     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     ethr =  5.22E-06,  avg # of iterations =  3.7
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96039  1.10303  6.06342
   spin  1
    eigenvalues: 
  0.987  0.987  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.179  0.388  0.000  0.140  0.293
  0.387  0.179  0.000  0.293  0.140
  0.405  0.028  0.000  0.532  0.035
  0.028  0.405  0.000  0.035  0.532
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.150  0.150  0.153  0.326  0.326
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.085  0.305  0.000  0.185  0.425
  0.305  0.085  0.000  0.425  0.185
  0.559  0.051  0.000  0.375  0.015
  0.051  0.559  0.000  0.015  0.375
    occupations:
  0.153  0.000  0.000  0.000  0.000
  0.000  0.257  0.000 -0.060  0.061
  0.000  0.000  0.257  0.061  0.060
  0.000 -0.060  0.061  0.218 -0.000
  0.000  0.061  0.060 -0.000  0.218
atomic mag. moment =     3.857358
atom    2   Tr[ns(na)] (up, down, total) =   1.16185  4.96478  6.12663
   spin  1
    eigenvalues: 
  0.161  0.161  0.165  0.338  0.338
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.362  0.031  0.000  0.605  0.002
  0.031  0.362  0.000  0.002  0.605
  0.466  0.142  0.000  0.220  0.172
  0.142  0.466  0.000  0.172  0.220
    occupations:
  0.165  0.000  0.000  0.000  0.000
  0.000  0.268  0.000  0.061  0.061
  0.000  0.000  0.268  0.061 -0.061
  0.000  0.061  0.061  0.230 -0.000
  0.000  0.061 -0.061 -0.000  0.230
   spin  2
    eigenvalues: 
  0.989  0.989  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.561  0.001  0.000  0.419  0.020
  0.001  0.561  0.000  0.020  0.419
  0.270  0.169  0.000  0.436  0.125
  0.169  0.270  0.000  0.125  0.436
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000 -0.003 -0.002
  0.000  0.000  0.992 -0.002  0.003
  0.000 -0.003 -0.002  0.993 -0.000
  0.000 -0.002  0.003 -0.000  0.993
atomic mag. moment =    -3.802934
atom    3   Tr[ns(na)] (up, down, total) =   4.96099  1.10821  6.06919
   spin  1
    eigenvalues: 
  0.987  0.987  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.182  0.393  0.000  0.132  0.293
  0.393  0.182  0.000  0.293  0.132
  0.399  0.025  0.000  0.540  0.036
  0.025  0.399  0.000  0.036  0.540
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.151  0.151  0.154  0.326  0.326
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.311  0.000  0.173  0.429
  0.311  0.087  0.000  0.429  0.173
  0.553  0.049  0.000  0.380  0.017
  0.049  0.553  0.000  0.017  0.380
    occupations:
  0.154  0.000  0.000  0.000  0.000
  0.000  0.256  0.000 -0.060  0.061
  0.000  0.000  0.256  0.061  0.060
  0.000 -0.060  0.061  0.221 -0.000
  0.000  0.061  0.060 -0.000  0.221
atomic mag. moment =     3.852781
atom    4   Tr[ns(na)] (up, down, total) =   1.12053  4.96250  6.08303
   spin  1
    eigenvalues: 
  0.153  0.153  0.156  0.329  0.329
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.362  0.031  0.000  0.604  0.003
  0.031  0.362  0.000  0.003  0.604
  0.465  0.141  0.000  0.221  0.173
  0.141  0.465  0.000  0.173  0.221
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.260  0.000  0.061  0.061
  0.000  0.000  0.260  0.061 -0.061
  0.000  0.061  0.061  0.222 -0.000
  0.000  0.061 -0.061 -0.000  0.222
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.569  0.001  0.000  0.407  0.023
  0.001  0.569  0.000  0.023  0.407
  0.261  0.169  0.000  0.444  0.126
  0.169  0.261  0.000  0.126  0.444
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.841972
atom    5   Tr[ns(na)] (up, down, total) =   4.96541  1.18138  6.14679
   spin  1
    eigenvalues: 
  0.989  0.989  0.996  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.189  0.393  0.000  0.128  0.290
  0.393  0.189  0.000  0.290  0.128
  0.395  0.023  0.000  0.545  0.038
  0.023  0.395  0.000  0.038  0.545
    occupations:
  0.996  0.000  0.000  0.000  0.000
  0.000  0.992  0.000  0.003 -0.003
  0.000  0.000  0.992 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.165  0.165  0.169  0.341  0.341
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.091  0.313  0.000  0.169  0.427
  0.313  0.091  0.000  0.427  0.169
  0.548  0.048  0.000  0.384  0.019
  0.048  0.548  0.000  0.019  0.384
    occupations:
  0.169  0.000  0.000  0.000  0.000
  0.000  0.270  0.000 -0.061  0.061
  0.000  0.000  0.270  0.061  0.061
  0.000 -0.061  0.061  0.236 -0.000
  0.000  0.061  0.061 -0.000  0.236
atomic mag. moment =     3.784027
atom    6   Tr[ns(na)] (up, down, total) =   1.11767  4.96031  6.07798
   spin  1
    eigenvalues: 
  0.153  0.153  0.156  0.328  0.328
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.372  0.033  0.000  0.592  0.003
  0.033  0.372  0.000  0.003  0.592
  0.454  0.141  0.000  0.228  0.176
  0.141  0.454  0.000  0.176  0.228
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.257  0.000  0.060  0.061
  0.000  0.000  0.257  0.061 -0.060
  0.000  0.060  0.061  0.224 -0.000
  0.000  0.061 -0.060 -0.000  0.224
   spin  2
    eigenvalues: 
  0.987  0.987  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.590  0.001  0.000  0.381  0.027
  0.001  0.590  0.000  0.027  0.381
  0.241  0.168  0.000  0.465  0.127
  0.168  0.241  0.000  0.127  0.465
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.990  0.000 -0.003 -0.003
  0.000  0.000  0.990 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.842634
N of occupied +U levels =   36.567035
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4649
 -------

     total cpu time spent up to now is     1491.0 secs

     total energy              =   -3013.09099630 Ry
     Harris-Foulkes estimate   =   -3013.09263358 Ry
     estimated scf accuracy    <       0.01051290 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    25.67 Bohr mag/cell

     iteration #  8     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     ethr =  3.58E-06,  avg # of iterations =  3.2
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96351  1.12169  6.08520
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.184  0.387  0.000  0.137  0.291
  0.387  0.184  0.000  0.291  0.137
  0.402  0.026  0.000  0.536  0.036
  0.026  0.402  0.000  0.036  0.536
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000  0.003 -0.003
  0.000  0.000  0.992 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.153  0.153  0.156  0.330  0.330
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.310  0.000  0.175  0.427
  0.310  0.088  0.000  0.427  0.175
  0.553  0.049  0.000  0.380  0.017
  0.049  0.553  0.000  0.017  0.380
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.260  0.000 -0.061  0.062
  0.000  0.000  0.260  0.062  0.061
  0.000 -0.061  0.062  0.223 -0.000
  0.000  0.062  0.061 -0.000  0.223
atomic mag. moment =     3.841828
atom    2   Tr[ns(na)] (up, down, total) =   1.11957  4.96335  6.08292
   spin  1
    eigenvalues: 
  0.153  0.153  0.156  0.329  0.329
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.032  0.000  0.600  0.003
  0.032  0.366  0.000  0.003  0.600
  0.461  0.141  0.000  0.224  0.174
  0.141  0.461  0.000  0.174  0.224
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.259  0.000  0.061  0.062
  0.000  0.000  0.259  0.062 -0.061
  0.000  0.061  0.062  0.223 -0.000
  0.000  0.062 -0.061 -0.000  0.223
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.404  0.022
  0.001  0.573  0.000  0.022  0.404
  0.258  0.168  0.000  0.448  0.126
  0.168  0.258  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.843786
atom    3   Tr[ns(na)] (up, down, total) =   4.96326  1.11887  6.08214
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.185  0.388  0.000  0.136  0.291
  0.388  0.185  0.000  0.291  0.136
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.152  0.152  0.155  0.329  0.329
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.176  0.427
  0.309  0.088  0.000  0.427  0.176
  0.553  0.049  0.000  0.380  0.017
  0.049  0.553  0.000  0.017  0.380
    occupations:
  0.155  0.000  0.000  0.000  0.000
  0.000  0.259  0.000 -0.061  0.062
  0.000  0.000  0.259  0.062  0.061
  0.000 -0.061  0.062  0.223 -0.000
  0.000  0.062  0.061 -0.000  0.223
atomic mag. moment =     3.844389
atom    4   Tr[ns(na)] (up, down, total) =   1.11424  4.96307  6.07731
   spin  1
    eigenvalues: 
  0.151  0.151  0.154  0.328  0.328
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.032  0.000  0.600  0.003
  0.032  0.365  0.000  0.003  0.600
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.154  0.000  0.000  0.000  0.000
  0.000  0.258  0.000  0.061  0.062
  0.000  0.000  0.258  0.062 -0.061
  0.000  0.061  0.062  0.222 -0.000
  0.000  0.062 -0.061 -0.000  0.222
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.572  0.001  0.000  0.405  0.022
  0.001  0.572  0.000  0.022  0.405
  0.260  0.168  0.000  0.446  0.126
  0.168  0.260  0.000  0.126  0.446
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.848828
atom    5   Tr[ns(na)] (up, down, total) =   4.96324  1.11779  6.08103
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.185  0.389  0.000  0.136  0.291
  0.389  0.185  0.000  0.291  0.136
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.152  0.152  0.155  0.329  0.329
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.309  0.000  0.176  0.427
  0.309  0.087  0.000  0.427  0.176
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.155  0.000  0.000  0.000  0.000
  0.000  0.259  0.000 -0.061  0.062
  0.000  0.000  0.259  0.062  0.061
  0.000 -0.061  0.062  0.222 -0.000
  0.000  0.062  0.061 -0.000  0.222
atomic mag. moment =     3.845448
atom    6   Tr[ns(na)] (up, down, total) =   1.12327  4.96358  6.08685
   spin  1
    eigenvalues: 
  0.153  0.153  0.156  0.330  0.330
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.031  0.000  0.601  0.003
  0.031  0.365  0.000  0.003  0.601
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.260  0.000  0.061  0.062
  0.000  0.000  0.260  0.062 -0.061
  0.000  0.061  0.062  0.223 -0.000
  0.000  0.062 -0.061 -0.000  0.223
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.571  0.001  0.000  0.406  0.022
  0.001  0.571  0.000  0.022  0.406
  0.260  0.168  0.000  0.446  0.126
  0.168  0.260  0.000  0.126  0.446
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000 -0.003 -0.003
  0.000  0.000  0.992 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.840314
N of occupied +U levels =   36.495447
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4666
 -------

     total cpu time spent up to now is     1670.7 secs

     total energy              =   -3013.09211175 Ry
     Harris-Foulkes estimate   =   -3013.09230971 Ry
     estimated scf accuracy    <       0.00087508 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.74 Bohr mag/cell

     iteration #  9     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     ethr =  2.98E-07,  avg # of iterations =  3.4
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96280  1.12927  6.09207
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.389  0.000  0.133  0.290
  0.389  0.187  0.000  0.290  0.133
  0.399  0.025  0.000  0.540  0.037
  0.025  0.399  0.000  0.037  0.540
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.332  0.332
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.310  0.000  0.174  0.428
  0.310  0.087  0.000  0.428  0.174
  0.553  0.049  0.000  0.380  0.017
  0.049  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.261  0.000 -0.061  0.062
  0.000  0.000  0.261  0.062  0.061
  0.000 -0.061  0.062  0.225 -0.000
  0.000  0.062  0.061 -0.000  0.225
atomic mag. moment =     3.833526
atom    2   Tr[ns(na)] (up, down, total) =   1.13019  4.96329  6.09348
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.333  0.333
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.031  0.000  0.601  0.003
  0.031  0.365  0.000  0.003  0.601
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.061  0.062
  0.000  0.000  0.262  0.062 -0.061
  0.000  0.061  0.062  0.225 -0.000
  0.000  0.062 -0.061 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.572  0.001  0.000  0.405  0.022
  0.001  0.572  0.000  0.022  0.405
  0.259  0.168  0.000  0.447  0.126
  0.168  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.833108
atom    3   Tr[ns(na)] (up, down, total) =   4.96357  1.13542  6.09899
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.185  0.386  0.000  0.138  0.291
  0.386  0.185  0.000  0.291  0.138
  0.402  0.026  0.000  0.536  0.036
  0.026  0.402  0.000  0.036  0.536
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.155  0.155  0.158  0.333  0.333
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.087  0.309  0.000  0.177  0.427
  0.309  0.087  0.000  0.427  0.177
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000 -0.061  0.062
  0.000  0.000  0.263  0.062  0.061
  0.000 -0.061  0.062  0.226 -0.000
  0.000  0.062  0.061 -0.000  0.226
atomic mag. moment =     3.828150
atom    4   Tr[ns(na)] (up, down, total) =   1.14369  4.96405  6.10774
   spin  1
    eigenvalues: 
  0.157  0.157  0.160  0.335  0.335
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.367  0.031  0.000  0.599  0.003
  0.031  0.367  0.000  0.003  0.599
  0.461  0.141  0.000  0.224  0.174
  0.141  0.461  0.000  0.174  0.224
    occupations:
  0.160  0.000  0.000  0.000  0.000
  0.000  0.264  0.000  0.061  0.062
  0.000  0.000  0.264  0.062 -0.061
  0.000  0.061  0.062  0.228 -0.000
  0.000  0.062 -0.061 -0.000  0.228
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.574  0.001  0.000  0.403  0.022
  0.001  0.574  0.000  0.022  0.403
  0.258  0.168  0.000  0.449  0.126
  0.168  0.258  0.000  0.126  0.449
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000 -0.003 -0.003
  0.000  0.000  0.992 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.820356
atom    5   Tr[ns(na)] (up, down, total) =   4.96374  1.13395  6.09769
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.186  0.387  0.000  0.136  0.290
  0.387  0.186  0.000  0.290  0.136
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.992  0.000  0.003 -0.003
  0.000  0.000  0.992 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.155  0.155  0.158  0.333  0.333
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.061  0.062
  0.000  0.000  0.262  0.062  0.061
  0.000 -0.061  0.062  0.226 -0.000
  0.000  0.062  0.061 -0.000  0.226
atomic mag. moment =     3.829781
atom    6   Tr[ns(na)] (up, down, total) =   1.12562  4.96284  6.08846
   spin  1
    eigenvalues: 
  0.153  0.153  0.156  0.331  0.331
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.032  0.000  0.601  0.003
  0.032  0.365  0.000  0.003  0.601
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.156  0.000  0.000  0.000  0.000
  0.000  0.261  0.000  0.061  0.062
  0.000  0.000  0.261  0.062 -0.061
  0.000  0.061  0.062  0.224 -0.000
  0.000  0.062 -0.061 -0.000  0.224
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.404  0.022
  0.001  0.573  0.000  0.022  0.404
  0.259  0.168  0.000  0.448  0.126
  0.168  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.837223
N of occupied +U levels =   36.578428
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4679
 -------

     total cpu time spent up to now is     1869.1 secs

     total energy              =   -3013.09234316 Ry
     Harris-Foulkes estimate   =   -3013.09237447 Ry
     estimated scf accuracy    <       0.00052259 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.70 Bohr mag/cell

     iteration # 10     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     ethr =  1.78E-07,  avg # of iterations =  3.1
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96334  1.13761  6.10095
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.136  0.290
  0.387  0.187  0.000  0.290  0.136
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.155  0.155  0.158  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000 -0.062  0.062
  0.000  0.000  0.263  0.062  0.062
  0.000 -0.062  0.062  0.226 -0.000
  0.000  0.062  0.062 -0.000  0.226
atomic mag. moment =     3.825729
atom    2   Tr[ns(na)] (up, down, total) =   1.13381  4.96319  6.09700
   spin  1
    eigenvalues: 
  0.155  0.155  0.158  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000  0.062  0.062
  0.000  0.000  0.263  0.062 -0.062
  0.000  0.062  0.062  0.226 -0.000
  0.000  0.062 -0.062 -0.000  0.226
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.574  0.001  0.000  0.403  0.022
  0.001  0.574  0.000  0.022  0.403
  0.258  0.167  0.000  0.449  0.126
  0.167  0.258  0.000  0.126  0.449
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.829371
atom    3   Tr[ns(na)] (up, down, total) =   4.96319  1.13531  6.09850
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.388  0.000  0.136  0.290
  0.388  0.187  0.000  0.290  0.136
  0.400  0.025  0.000  0.538  0.036
  0.025  0.400  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.155  0.155  0.158  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.176  0.427
  0.309  0.088  0.000  0.427  0.176
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000 -0.062  0.062
  0.000  0.000  0.263  0.062  0.062
  0.000 -0.062  0.062  0.226 -0.000
  0.000  0.062  0.062 -0.000  0.226
atomic mag. moment =     3.827876
atom    4   Tr[ns(na)] (up, down, total) =   1.13365  4.96308  6.09673
   spin  1
    eigenvalues: 
  0.154  0.154  0.158  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.031  0.000  0.600  0.003
  0.031  0.365  0.000  0.003  0.600
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000  0.062  0.062
  0.000  0.000  0.263  0.062 -0.062
  0.000  0.062  0.062  0.226 -0.000
  0.000  0.062 -0.062 -0.000  0.226
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.404  0.022
  0.001  0.573  0.000  0.022  0.404
  0.258  0.168  0.000  0.448  0.126
  0.168  0.258  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.829433
atom    5   Tr[ns(na)] (up, down, total) =   4.96311  1.13144  6.09455
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.136  0.290
  0.387  0.187  0.000  0.290  0.136
  0.401  0.026  0.000  0.538  0.036
  0.026  0.401  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.333  0.333
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.427
  0.309  0.088  0.000  0.427  0.177
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.831672
atom    6   Tr[ns(na)] (up, down, total) =   1.13617  4.96332  6.09949
   spin  1
    eigenvalues: 
  0.155  0.155  0.158  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.223  0.174
  0.141  0.462  0.000  0.174  0.223
    occupations:
  0.158  0.000  0.000  0.000  0.000
  0.000  0.263  0.000  0.062  0.062
  0.000  0.000  0.263  0.062 -0.062
  0.000  0.062  0.062  0.226 -0.000
  0.000  0.062 -0.062 -0.000  0.226
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.168  0.000  0.447  0.126
  0.168  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.827142
N of occupied +U levels =   36.587219
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4695
 -------

     total cpu time spent up to now is     2042.4 secs

     total energy              =   -3013.09239474 Ry
     Harris-Foulkes estimate   =   -3013.09240275 Ry
     estimated scf accuracy    <       0.00019046 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.69 Bohr mag/cell

     iteration # 11     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     ethr =  6.48E-08,  avg # of iterations =  3.1
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96284  1.13227  6.09511
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.136  0.290
  0.387  0.187  0.000  0.290  0.136
  0.400  0.026  0.000  0.538  0.036
  0.026  0.400  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.049  0.000  0.380  0.017
  0.049  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.830577
atom    2   Tr[ns(na)] (up, down, total) =   1.13253  4.96290  6.09542
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.461  0.141  0.000  0.224  0.174
  0.141  0.461  0.000  0.174  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.062
  0.000  0.000  0.262  0.062 -0.062
  0.000  0.062  0.062  0.225 -0.000
  0.000  0.062 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.404  0.022
  0.001  0.573  0.000  0.022  0.404
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.830371
atom    3   Tr[ns(na)] (up, down, total) =   4.96294  1.13303  6.09596
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.137  0.290
  0.387  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.178  0.426
  0.309  0.088  0.000  0.426  0.178
  0.554  0.050  0.000  0.380  0.017
  0.050  0.554  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.829906
atom    4   Tr[ns(na)] (up, down, total) =   1.13332  4.96300  6.09632
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.174
  0.141  0.462  0.000  0.174  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.062
  0.000  0.000  0.262  0.062 -0.062
  0.000  0.062  0.062  0.226 -0.000
  0.000  0.062 -0.062 -0.000  0.226
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.404  0.022
  0.001  0.573  0.000  0.022  0.404
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.829681
atom    5   Tr[ns(na)] (up, down, total) =   4.96293  1.13276  6.09569
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.136  0.290
  0.387  0.187  0.000  0.290  0.136
  0.401  0.026  0.000  0.538  0.036
  0.026  0.401  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.830163
atom    6   Tr[ns(na)] (up, down, total) =   1.13288  4.96288  6.09576
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.365  0.031  0.000  0.601  0.003
  0.031  0.365  0.000  0.003  0.601
  0.462  0.141  0.000  0.223  0.173
  0.141  0.462  0.000  0.173  0.223
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.263  0.000  0.062  0.062
  0.000  0.000  0.263  0.062 -0.062
  0.000  0.062  0.062  0.225 -0.000
  0.000  0.062 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.830007
N of occupied +U levels =   36.574264
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4697
 -------

     total cpu time spent up to now is     2218.1 secs

     total energy              =   -3013.09241223 Ry
     Harris-Foulkes estimate   =   -3013.09241432 Ry
     estimated scf accuracy    <       0.00002334 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    25.70 Bohr mag/cell

     iteration # 12     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     c_bands:  1 eigenvalues not converged
     ethr =  7.94E-09,  avg # of iterations =  3.4
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96291  1.13238  6.09529
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.137  0.290
  0.387  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.830532
atom    2   Tr[ns(na)] (up, down, total) =   1.13126  4.96286  6.09412
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.223  0.173
  0.141  0.462  0.000  0.173  0.223
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.062
  0.000  0.000  0.262  0.062 -0.062
  0.000  0.062  0.062  0.225 -0.000
  0.000  0.062 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.447  0.126
  0.167  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.831603
atom    3   Tr[ns(na)] (up, down, total) =   4.96290  1.13208  6.09498
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.137  0.290
  0.387  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.538  0.036
  0.026  0.401  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.830824
atom    4   Tr[ns(na)] (up, down, total) =   1.13149  4.96284  6.09432
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.174
  0.141  0.462  0.000  0.174  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.831346
atom    5   Tr[ns(na)] (up, down, total) =   4.96291  1.13139  6.09430
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.062
  0.000  0.000  0.262  0.062  0.062
  0.000 -0.062  0.062  0.225 -0.000
  0.000  0.062  0.062 -0.000  0.225
atomic mag. moment =     3.831522
atom    6   Tr[ns(na)] (up, down, total) =   1.13172  4.96288  6.09461
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.062
  0.000  0.000  0.262  0.062 -0.062
  0.000  0.062  0.062  0.225 -0.000
  0.000  0.062 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.831162
N of occupied +U levels =   36.567626
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4694
 -------

     total cpu time spent up to now is     2416.6 secs

     total energy              =   -3013.09241804 Ry
     Harris-Foulkes estimate   =   -3013.09241855 Ry
     estimated scf accuracy    <       0.00000790 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.70 Bohr mag/cell

     iteration # 13     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  2.69E-09,  avg # of iterations =  3.3
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96288  1.13183  6.09471
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831048
atom    2   Tr[ns(na)] (up, down, total) =   1.13248  4.96296  6.09544
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.830485
atom    3   Tr[ns(na)] (up, down, total) =   4.96288  1.13146  6.09434
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831416
atom    4   Tr[ns(na)] (up, down, total) =   1.13253  4.96298  6.09551
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.447  0.126
  0.167  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.830447
atom    5   Tr[ns(na)] (up, down, total) =   4.96289  1.13153  6.09442
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.387  0.000  0.137  0.290
  0.387  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.538  0.036
  0.026  0.401  0.000  0.036  0.538
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831365
atom    6   Tr[ns(na)] (up, down, total) =   1.13173  4.96290  6.09463
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.831175
N of occupied +U levels =   36.569057
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4694
 -------

     total cpu time spent up to now is     2610.0 secs

     total energy              =   -3013.09241966 Ry
     Harris-Foulkes estimate   =   -3013.09241990 Ry
     estimated scf accuracy    <       0.00000340 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    25.70 Bohr mag/cell

     iteration # 14     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  1.16E-09,  avg # of iterations =  3.2
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96290  1.13187  6.09477
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831033
atom    2   Tr[ns(na)] (up, down, total) =   1.13149  4.96286  6.09435
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.447  0.126
  0.167  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.831371
atom    3   Tr[ns(na)] (up, down, total) =   4.96294  1.13232  6.09526
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.830617
atom    4   Tr[ns(na)] (up, down, total) =   1.13149  4.96287  6.09436
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.461  0.141  0.000  0.224  0.173
  0.141  0.461  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.831371
atom    5   Tr[ns(na)] (up, down, total) =   4.96296  1.13247  6.09543
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.830484
atom    6   Tr[ns(na)] (up, down, total) =   1.13197  4.96288  6.09486
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.992 -0.000
  0.000 -0.003  0.003 -0.000  0.992
atomic mag. moment =    -3.830910
N of occupied +U levels =   36.569018
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4693
 -------

     total cpu time spent up to now is     2792.0 secs

     total energy              =   -3013.09242044 Ry
     Harris-Foulkes estimate   =   -3013.09242049 Ry
     estimated scf accuracy    <       0.00000191 Ry

     total magnetization       =    -0.00 Bohr mag/cell
     absolute magnetization    =    25.70 Bohr mag/cell

     iteration # 15     ecut=    40.00 Ry     beta= 0.40
     CG style diagonalization
     ethr =  6.48E-10,  avg # of iterations =  3.1
 --- enter write_ns ---
 LDA+U parameters:
U( 3)     =  2.43461000
alpha( 3) =  0.00000000
U( 4)     =  2.43461000
alpha( 4) =  0.00000000
atom    1   Tr[ns(na)] (up, down, total) =   4.96289  1.13178  6.09467
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.992 -0.000
  0.000 -0.003 -0.003 -0.000  0.992
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831110
atom    2   Tr[ns(na)] (up, down, total) =   1.13195  4.96291  6.09486
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.447  0.126
  0.167  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.830968
atom    3   Tr[ns(na)] (up, down, total) =   4.96289  1.13172  6.09462
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831174
atom    4   Tr[ns(na)] (up, down, total) =   1.13198  4.96291  6.09489
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.461  0.141  0.000  0.224  0.173
  0.141  0.461  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.447  0.126
  0.167  0.259  0.000  0.126  0.447
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.830935
atom    5   Tr[ns(na)] (up, down, total) =   4.96291  1.13172  6.09463
   spin  1
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.187  0.386  0.000  0.137  0.290
  0.386  0.187  0.000  0.290  0.137
  0.401  0.026  0.000  0.537  0.036
  0.026  0.401  0.000  0.036  0.537
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000  0.003 -0.003
  0.000  0.000  0.991 -0.003 -0.003
  0.000  0.003 -0.003  0.993 -0.000
  0.000 -0.003 -0.003 -0.000  0.993
   spin  2
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.088  0.309  0.000  0.177  0.426
  0.309  0.088  0.000  0.426  0.177
  0.553  0.050  0.000  0.380  0.017
  0.050  0.553  0.000  0.017  0.380
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000 -0.062  0.063
  0.000  0.000  0.262  0.063  0.062
  0.000 -0.062  0.063  0.225 -0.000
  0.000  0.063  0.062 -0.000  0.225
atomic mag. moment =     3.831182
atom    6   Tr[ns(na)] (up, down, total) =   1.13190  4.96291  6.09481
   spin  1
    eigenvalues: 
  0.154  0.154  0.157  0.334  0.334
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.366  0.031  0.000  0.600  0.003
  0.031  0.366  0.000  0.003  0.600
  0.462  0.141  0.000  0.224  0.173
  0.141  0.462  0.000  0.173  0.224
    occupations:
  0.157  0.000  0.000  0.000  0.000
  0.000  0.262  0.000  0.062  0.063
  0.000  0.000  0.262  0.063 -0.062
  0.000  0.062  0.063  0.225 -0.000
  0.000  0.063 -0.062 -0.000  0.225
   spin  2
    eigenvalues: 
  0.988  0.988  0.995  0.996  0.996
    eigenvectors:
  0.000  0.000  1.000  0.000  0.000
  0.573  0.001  0.000  0.405  0.022
  0.001  0.573  0.000  0.022  0.405
  0.259  0.167  0.000  0.448  0.126
  0.167  0.259  0.000  0.126  0.448
    occupations:
  0.995  0.000  0.000  0.000  0.000
  0.000  0.991  0.000 -0.003 -0.003
  0.000  0.000  0.991 -0.003  0.003
  0.000 -0.003 -0.003  0.993 -0.000
  0.000 -0.003  0.003 -0.000  0.993
atomic mag. moment =    -3.831012
N of occupied +U levels =   36.568467
 --- exit write_ns ---
 --- in v_hubbard ---
Hubbard energy    0.4691
 -------

     Magnetic moment per site:
     atom:    1    charge:   12.9983    magn:    3.7010    constr:    0.0000
     atom:    2    charge:   12.9989    magn:   -3.7011    constr:    0.0000
     atom:    3    charge:   12.9986    magn:    3.7012    constr:    0.0000
     atom:    4    charge:   12.9991    magn:   -3.7011    constr:    0.0000
     atom:    5    charge:   12.9989    magn:    3.7013    constr:    0.0000
     atom:    6    charge:   12.9984    magn:   -3.7010    constr:    0.0000
     atom:    7    charge:   10.5985    magn:   -0.0005    constr:    0.0000
     atom:    8    charge:   10.5986    magn:    0.0005    constr:    0.0000
     atom:    9    charge:   10.5986    magn:   -0.0005    constr:    0.0000
     atom:   10    charge:   10.5988    magn:    0.0005    constr:    0.0000
     atom:   11    charge:   10.5988    magn:   -0.0005    constr:    0.0000
     atom:   12    charge:   10.5985    magn:    0.0005    constr:    0.0000
     atom:   13    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   14    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   15    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   16    charge:    5.3694    magn:   -0.0609    constr:    0.0000
     atom:   17    charge:    5.3694    magn:   -0.0609    constr:    0.0000
     atom:   18    charge:    5.3694    magn:   -0.0609    constr:    0.0000
     atom:   19    charge:    5.3693    magn:    0.0609    constr:    0.0000
     atom:   20    charge:    5.3693    magn:    0.0609    constr:    0.0000
     atom:   21    charge:    5.3693    magn:    0.0609    constr:    0.0000
     atom:   22    charge:    5.3691    magn:   -0.0609    constr:    0.0000
     atom:   23    charge:    5.3691    magn:   -0.0609    constr:    0.0000
     atom:   24    charge:    5.3691    magn:   -0.0609    constr:    0.0000
     atom:   25    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   26    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   27    charge:    5.3692    magn:    0.0609    constr:    0.0000
     atom:   28    charge:    5.3692    magn:   -0.0609    constr:    0.0000
     atom:   29    charge:    5.3692    magn:   -0.0609    constr:    0.0000
     atom:   30    charge:    5.3692    magn:   -0.0609    constr:    0.0000

     total cpu time spent up to now is     2963.5 secs

     End of self-consistent calculation

 ------ SPIN UP ------------


          k = 0.0000 0.0000 0.0000 ( 11177 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1990 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0367 -12.0367 -12.0015 -12.0015 -12.0011 -12.0011 -11.9479 -11.9473
   -11.9473 -11.9471 -11.9470 -11.9446 -11.9446 -11.9349 -11.9348 -11.9049
   -11.9049 -11.9047 -11.9047 -11.9030 -11.9028 -11.9008 -11.9008 -11.8893
   -11.8792 -11.8792 -11.8791 -11.8791 -11.8571 -11.8571  -8.2663  -7.9016
    -7.9016  -7.0598  -7.0598  -6.5738  -6.4976  -6.4976  -6.4733  -6.4733
    -6.4425  -6.4423  -6.2652  -6.2651  -6.2325  -6.2324  -6.1709  -6.1709
    -0.0356   0.1704   0.1704   0.9409   0.9410   1.6475   3.8944   3.8944
     4.1522   4.1522   4.1630   4.1681   4.1681   4.4056   4.4056   4.4354
     4.4354   4.4750   4.4755   4.5661   4.5661   4.6165   5.0433   5.0435
     5.4242   5.4242   5.4499   5.5974   5.5974   5.6319   5.6320   6.0637
     6.0639   6.3688   6.3688   6.4092   6.4092   6.5986   6.8182   6.8182
     7.0755   7.0756   7.3686   7.3686   7.7263   7.7264   7.7428   7.7428
     7.7977   7.7977   7.8772   7.8774   7.9407   7.9408   8.5477   8.5477
     8.5536   8.5724   8.5724   8.6478   8.6478   8.9700   8.9700   9.1433
     9.3406   9.3410   9.4044   9.4044   9.4368   9.6827   9.6827   9.7339
     9.7344   9.8690   9.8690  11.7185  11.7189  11.7504  11.7504  11.7801
    11.7836  11.7836  11.8593  11.8593  12.9428  12.9428  13.0448  13.0449
    13.2263  13.2263  13.9831  13.9835  14.1286  14.1291  15.0461  15.0462
    15.1435  15.1437  15.1657  15.1658  15.5720  15.7169  15.7206  15.8614

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.0000 0.1335 ( 11196 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1990 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0327 -12.0324 -12.0193 -12.0190 -11.9819 -11.9815 -11.9612 -11.9608
   -11.9478 -11.9440 -11.9439 -11.9425 -11.9422 -11.9392 -11.9227 -11.9135
   -11.9039 -11.9038 -11.9027 -11.9026 -11.9002 -11.9000 -11.8941 -11.8908
   -11.8906 -11.8890 -11.8692 -11.8692 -11.8626 -11.8625  -8.2237  -8.0990
    -7.6512  -7.3447  -6.8085  -6.6357  -6.5004  -6.4923  -6.4851  -6.4650
    -6.4480  -6.4322  -6.2859  -6.2695  -6.2255  -6.1962  -6.1900  -6.1734
    -0.0161   0.0483   0.3602   0.6260   1.2713   1.5402   3.9389   3.9409
     4.0503   4.0592   4.1940   4.1965   4.2042   4.3231   4.3464   4.3786
     4.4810   4.5131   4.5492   4.5769   4.6130   4.6759   4.8326   5.2439
     5.4272   5.4955   5.5069   5.5729   5.5825   5.6345   5.6401   5.8484
     6.0662   6.3428   6.3627   6.4938   6.5067   6.6360   6.6681   6.7746
     6.8467   7.2182   7.3356   7.4426   7.6831   7.6972   7.7066   7.7264
     7.7668   7.9397   7.9908   8.0377   8.1058   8.1357   8.3536   8.3605
     8.4158   8.4187   8.4234   8.7802   8.8284   8.9062   8.9347   9.1760
     9.2525   9.4060   9.4079   9.4888   9.5304   9.5813   9.5978   9.7631
     9.8150   9.8565   9.8852  11.7024  11.7315  11.7451  11.7468  11.7699
    11.7951  11.8236  11.8388  11.8523  12.9047  12.9458  13.0425  13.1547
    13.1553  13.2261  13.5217  13.7108  14.6688  14.8046  15.0363  15.0580
    15.0618  15.0913  15.2917  15.2978  15.3013  15.5044  15.5091  16.0725

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.3849 0.0000 ( 11182 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0278 -12.0191 -12.0038 -11.9969 -11.9816 -11.9789 -11.9651 -11.9589
   -11.9585 -11.9422 -11.9363 -11.9352 -11.9315 -11.9262 -11.9180 -11.9159
   -11.9103 -11.9067 -11.9062 -11.9057 -11.9048 -11.8994 -11.8960 -11.8932
   -11.8923 -11.8891 -11.8824 -11.8778 -11.8703 -11.8648  -7.9298  -7.6611
    -7.6125  -7.2183  -6.9747  -6.9668  -6.6881  -6.6055  -6.5898  -6.5510
    -6.4995  -6.3826  -6.3585  -6.3142  -6.2991  -6.2492  -6.2092  -6.1886
     0.0844   0.1032   0.3168   0.5798   0.8072   0.9918   3.8684   4.0413
     4.0753   4.2194   4.2413   4.2965   4.3696   4.4450   4.4572   4.4680
     4.5984   4.6162   4.6862   4.8047   4.8755   4.9989   5.0441   5.0940
     5.3674   5.6678   5.7681   5.7755   5.8993   5.9813   6.0052   6.3136
     6.3574   6.3899   6.4866   6.5434   6.6929   6.7893   6.8224   6.9036
     7.0890   7.3218   7.3511   7.3854   7.3935   7.4864   7.5095   7.6795
     7.7014   7.7789   7.9328   7.9926   8.0790   8.2021   8.2449   8.3066
     8.3295   8.4029   8.4062   8.4968   8.5328   8.6040   8.6821   8.8599
     8.9853   9.0891   9.3776   9.3870   9.4620   9.5350   9.6819   9.7596
     9.8796   9.8880   9.9521  11.7435  11.7617  11.7623  11.7971  11.8340
    11.8762  11.9406  11.9514  11.9777  12.4924  12.7699  12.8660  12.9059
    13.0422  13.1524  13.4714  13.9559  14.0039  14.0574  14.3466  14.4101
    14.4843  14.5149  14.9640  15.3032  15.3783  15.3914  15.4665  15.5452

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.3849 0.1335 ( 11170 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0284 -12.0168 -12.0067 -11.9880 -11.9853 -11.9852 -11.9640 -11.9631
   -11.9486 -11.9475 -11.9353 -11.9339 -11.9291 -11.9273 -11.9206 -11.9130
   -11.9095 -11.9092 -11.9070 -11.9041 -11.9023 -11.9003 -11.8997 -11.8919
   -11.8885 -11.8871 -11.8866 -11.8757 -11.8708 -11.8631  -7.9003  -7.7790
    -7.4963  -7.1646  -7.1370  -6.8912  -6.7196  -6.6045  -6.5925  -6.4930
    -6.4522  -6.4270  -6.3486  -6.3376  -6.3331  -6.2109  -6.2057  -6.2053
     0.0511   0.2088   0.2188   0.6375   0.7635   1.0035   3.8923   3.9575
     4.1673   4.2022   4.2214   4.3469   4.3486   4.3759   4.4719   4.4823
     4.5786   4.6394   4.8006   4.8059   4.8448   4.9795   4.9855   5.1098
     5.3878   5.5805   5.6946   5.7758   5.8492   6.0738   6.1673   6.1851
     6.4119   6.4160   6.4765   6.5494   6.7215   6.7712   6.7968   6.9666
     7.1017   7.2379   7.2933   7.3725   7.4225   7.4942   7.5384   7.7065
     7.7697   7.8047   7.8711   7.8731   8.1367   8.2005   8.2126   8.3359
     8.3601   8.4123   8.4411   8.4624   8.4970   8.5684   8.6742   8.9356
     8.9557   9.1184   9.2690   9.4198   9.5290   9.5586   9.6538   9.7478
     9.8374   9.9390   9.9502  11.7349  11.7699  11.7869  11.7876  11.8110
    11.8851  11.9414  11.9554  11.9668  12.5311  12.7133  12.8166  12.9610
    13.0391  13.1570  13.5219  13.8738  13.9426  14.1196  14.1769  14.4113
    14.6623  14.7687  14.8018  15.2313  15.3140  15.4517  15.5022  15.5841

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333 0.5774 0.0000 ( 11157 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0015 -12.0014 -12.0014 -11.9977 -11.9977 -11.9977 -11.9520 -11.9519
   -11.9519 -11.9411 -11.9410 -11.9409 -11.9178 -11.9177 -11.9177 -11.9128
   -11.9128 -11.9127 -11.9034 -11.9033 -11.9033 -11.9000 -11.8999 -11.8999
   -11.8847 -11.8847 -11.8847 -11.8800 -11.8800 -11.8799  -7.4098  -7.4092
    -7.4091  -7.3147  -7.3146  -7.3140  -6.7071  -6.7070  -6.7068  -6.4854
    -6.4853  -6.4852  -6.3680  -6.3679  -6.3677  -6.2136  -6.2134  -6.2132
     0.3867   0.3868   0.3869   0.4415   0.4415   0.4416   3.9058   3.9065
     3.9066   4.3077   4.3081   4.3081   4.4774   4.4778   4.4778   4.8352
     4.8354   4.8356   4.9244   4.9250   4.9253   5.0556   5.0556   5.0561
     5.6954   5.6954   5.6957   5.9339   5.9339   5.9341   6.2650   6.2652
     6.2654   6.4699   6.4699   6.4700   6.8227   6.8227   6.8228   7.0552
     7.0552   7.0554   7.3279   7.3281   7.3281   7.6287   7.6288   7.6290
     7.7872   7.7876   7.7876   7.9716   7.9718   7.9720   8.3568   8.3571
     8.3571   8.5086   8.5086   8.5087   8.6103   8.6104   8.6105   9.1419
     9.1421   9.1423   9.2826   9.2828   9.2830   9.4649   9.4652   9.4654
     9.7928   9.7932   9.7935  11.8072  11.8074  11.8075  11.9077  11.9078
    11.9078  11.9845  11.9847  11.9850  12.6073  12.6074  12.6076  12.9886
    12.9889  12.9889  13.6158  13.6159  13.6160  13.8818  13.8819  13.8820
    14.7850  14.7850  14.7852  14.9893  14.9893  14.9898  15.4346  15.4349

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333 0.5774 0.1335 ( 11151 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0007 -12.0006 -12.0005 -11.9983 -11.9982 -11.9982 -11.9497 -11.9496
   -11.9496 -11.9438 -11.9437 -11.9436 -11.9160 -11.9159 -11.9158 -11.9133
   -11.9133 -11.9132 -11.9030 -11.9029 -11.9029 -11.9010 -11.9009 -11.9009
   -11.8831 -11.8831 -11.8831 -11.8809 -11.8809 -11.8808  -7.4033  -7.4025
    -7.4024  -7.3245  -7.3245  -7.3237  -6.6719  -6.6718  -6.6716  -6.5523
    -6.5523  -6.5521  -6.3146  -6.3144  -6.3141  -6.2322  -6.2320  -6.2317
     0.4002   0.4003   0.4003   0.4276   0.4276   0.4278   3.9004   3.9011
     3.9011   4.3563   4.3567   4.3567   4.4451   4.4454   4.4455   4.8203
     4.8204   4.8204   4.9293   4.9297   4.9300   5.0527   5.0527   5.0532
     5.6807   5.6808   5.6810   5.9292   5.9293   5.9295   6.2631   6.2632
     6.2633   6.5158   6.5158   6.5159   6.8218   6.8219   6.8219   7.1133
     7.1135   7.1136   7.2236   7.2238   7.2239   7.6052   7.6053   7.6055
     7.8240   7.8240   7.8241   8.0024   8.0025   8.0026   8.3359   8.3361
     8.3361   8.4826   8.4826   8.4827   8.6485   8.6485   8.6486   9.1593
     9.1595   9.1598   9.2571   9.2574   9.2576   9.4514   9.4516   9.4518
     9.8020   9.8023   9.8027  11.8080  11.8082  11.8083  11.9141  11.9142
    11.9143  11.9808  11.9810  11.9813  12.6329  12.6331  12.6333  12.9669
    12.9672  12.9673  13.6223  13.6225  13.6225  13.8225  13.8226  13.8228
    14.8794  14.8798  14.8802  14.9908  14.9910  14.9912  15.4817  15.4836

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k =-0.3333 0.1925 0.1335 ( 11170 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0311 -12.0047 -12.0033 -12.0006 -11.9968 -11.9720 -11.9692 -11.9542
   -11.9496 -11.9489 -11.9382 -11.9366 -11.9259 -11.9242 -11.9218 -11.9140
   -11.9105 -11.9097 -11.9061 -11.9060 -11.9041 -11.9005 -11.8961 -11.8929
   -11.8909 -11.8835 -11.8805 -11.8791 -11.8773 -11.8606  -7.8885  -7.8056
    -7.4247  -7.3065  -7.0546  -6.8746  -6.6973  -6.6474  -6.5553  -6.5255
    -6.4964  -6.3918  -6.3790  -6.3126  -6.2957  -6.2318  -6.2221  -6.1889
     0.0502   0.1356   0.3853   0.4657   0.9078   0.9385   3.8952   3.9669
     4.1338   4.2000   4.2885   4.3013   4.3456   4.4156   4.4383   4.5241
     4.5327   4.6655   4.7106   4.8052   4.8800   4.9794   4.9799   5.1257
     5.4708   5.5577   5.6846   5.7503   5.8475   6.0394   6.1777   6.2978
     6.2980   6.3897   6.4801   6.5506   6.7843   6.7989   6.8188   6.8694
     7.1537   7.1815   7.3620   7.3939   7.4586   7.4889   7.5752   7.5802
     7.7557   7.7741   7.7929   8.0341   8.0977   8.1779   8.2315   8.3339
     8.3468   8.3845   8.4784   8.4910   8.5364   8.6039   8.6580   8.8378
     8.9191   9.1319   9.2764   9.4689   9.5066   9.5742   9.6899   9.6954
     9.8810   9.9026   9.9519  11.7456  11.7502  11.7704  11.7903  11.8501
    11.8768  11.9327  11.9569  11.9706  12.5637  12.6539  12.8992  12.9411
    12.9699  13.2069  13.6239  13.7485  13.9708  14.1050  14.2691  14.3243
    14.5932  14.6421  15.1567  15.2359  15.2388  15.2795  15.3278  15.6809

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333-0.5774-0.1335 ( 11151 PWs)   bands (ev):

   -79.4733 -79.4726 -79.4722 -74.6889 -74.6887 -74.6884 -45.8700 -45.8693
   -45.8690 -45.8099 -45.8099 -45.8092 -45.8092 -45.8088 -45.8088 -41.2483
   -41.2481 -41.2478 -41.1991 -41.1990 -41.1988 -41.1988 -41.1985 -41.1985
   -12.0005 -12.0004 -12.0003 -11.9984 -11.9984 -11.9983 -11.9495 -11.9494
   -11.9493 -11.9441 -11.9439 -11.9438 -11.9157 -11.9155 -11.9155 -11.9136
   -11.9135 -11.9134 -11.9029 -11.9028 -11.9027 -11.9011 -11.9011 -11.9010
   -11.8832 -11.8831 -11.8831 -11.8808 -11.8808 -11.8807  -7.4039  -7.4032
    -7.4031  -7.3236  -7.3235  -7.3228  -6.6768  -6.6766  -6.6765  -6.5449
    -6.5448  -6.5446  -6.3201  -6.3199  -6.3197  -6.2295  -6.2293  -6.2290
     0.3998   0.3999   0.3999   0.4281   0.4281   0.4282   3.8962   3.8968
     3.8969   4.3594   4.3598   4.3599   4.4491   4.4494   4.4495   4.8170
     4.8171   4.8172   4.9353   4.9357   4.9360   5.0567   5.0568   5.0573
     5.6391   5.6392   5.6394   5.9624   5.9626   5.9627   6.2780   6.2781
     6.2782   6.5018   6.5019   6.5020   6.8350   6.8351   6.8352   7.0639
     7.0640   7.0643   7.2468   7.2469   7.2470   7.6329   7.6331   7.6331
     7.7584   7.7585   7.7585   8.0703   8.0705   8.0705   8.3701   8.3704
     8.3704   8.4812   8.4812   8.4813   8.5892   8.5892   8.5894   9.1543
     9.1545   9.1547   9.2467   9.2469   9.2471   9.4722   9.4725   9.4726
     9.8063   9.8066   9.8070  11.8085  11.8087  11.8088  11.8981  11.8981
    11.8982  11.9919  11.9920  11.9923  12.6104  12.6106  12.6108  12.9948
    12.9951  12.9951  13.6732  13.6735  13.6736  13.7680  13.7681  13.7684
    14.8818  14.8819  14.8820  14.9807  14.9807  14.9809  15.4994  15.4995

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

 ------ SPIN DOWN ----------


          k = 0.0000 0.0000 0.0000 ( 11177 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8100 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0367 -12.0367 -12.0015 -12.0015 -12.0011 -12.0011 -11.9479 -11.9473
   -11.9473 -11.9471 -11.9471 -11.9446 -11.9446 -11.9349 -11.9348 -11.9048
   -11.9048 -11.9047 -11.9047 -11.9029 -11.9029 -11.9008 -11.9008 -11.8893
   -11.8792 -11.8792 -11.8791 -11.8791 -11.8571 -11.8571  -8.2663  -7.9016
    -7.9015  -7.0599  -7.0597  -6.5738  -6.4976  -6.4976  -6.4732  -6.4732
    -6.4424  -6.4422  -6.2653  -6.2652  -6.2326  -6.2325  -6.1709  -6.1709
    -0.0356   0.1704   0.1704   0.9409   0.9410   1.6475   3.8942   3.8942
     4.1519   4.1519   4.1630   4.1678   4.1678   4.4053   4.4053   4.4351
     4.4351   4.4750   4.4752   4.5657   4.5657   4.6161   5.0431   5.0434
     5.4241   5.4241   5.4499   5.5974   5.5974   5.6319   5.6319   6.0635
     6.0639   6.3686   6.3686   6.4090   6.4090   6.5984   6.8181   6.8182
     7.0754   7.0756   7.3685   7.3685   7.7264   7.7265   7.7428   7.7428
     7.7977   7.7977   7.8772   7.8774   7.9407   7.9408   8.5477   8.5477
     8.5536   8.5724   8.5724   8.6477   8.6477   8.9701   8.9701   9.1434
     9.3407   9.3412   9.4044   9.4044   9.4367   9.6826   9.6826   9.7338
     9.7344   9.8689   9.8690  11.7188  11.7193  11.7508  11.7508  11.7805
    11.7840  11.7840  11.8597  11.8597  12.9431  12.9431  13.0451  13.0452
    13.2265  13.2266  13.9830  13.9835  14.1287  14.1291  15.0461  15.0462
    15.1435  15.1440  15.1657  15.1659  15.5683  15.7162  15.7377  15.8620

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.0000 0.1335 ( 11196 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8100 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0327 -12.0324 -12.0193 -12.0190 -11.9819 -11.9815 -11.9612 -11.9608
   -11.9478 -11.9440 -11.9439 -11.9425 -11.9422 -11.9392 -11.9227 -11.9135
   -11.9039 -11.9038 -11.9027 -11.9026 -11.9002 -11.9000 -11.8941 -11.8908
   -11.8906 -11.8890 -11.8692 -11.8692 -11.8626 -11.8625  -8.2237  -8.0990
    -7.6512  -7.3447  -6.8086  -6.6357  -6.5004  -6.4922  -6.4851  -6.4649
    -6.4479  -6.4321  -6.2861  -6.2696  -6.2256  -6.1963  -6.1900  -6.1734
    -0.0161   0.0483   0.3602   0.6261   1.2714   1.5402   3.9387   3.9406
     4.0500   4.0590   4.1937   4.1962   4.2042   4.3230   4.3461   4.3783
     4.4806   4.5127   4.5488   4.5766   4.6126   4.6755   4.8323   5.2439
     5.4271   5.4955   5.5069   5.5728   5.5824   5.6345   5.6401   5.8484
     6.0659   6.3426   6.3624   6.4936   6.5066   6.6359   6.6679   6.7746
     6.8467   7.2181   7.3356   7.4426   7.6831   7.6973   7.7067   7.7264
     7.7668   7.9398   7.9908   8.0377   8.1059   8.1357   8.3535   8.3605
     8.4158   8.4186   8.4233   8.7803   8.8284   8.9062   8.9348   9.1761
     9.2526   9.4061   9.4079   9.4887   9.5304   9.5813   9.5977   9.7631
     9.8149   9.8565   9.8851  11.7027  11.7319  11.7454  11.7472  11.7702
    11.7955  11.8239  11.8391  11.8526  12.9050  12.9461  13.0427  13.1548
    13.1555  13.2264  13.5216  13.7109  14.6687  14.8046  15.0362  15.0578
    15.0613  15.0913  15.2917  15.2975  15.3000  15.5001  15.5091  15.5151

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.3849 0.0000 ( 11182 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8101 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0278 -12.0191 -12.0038 -11.9969 -11.9816 -11.9789 -11.9651 -11.9589
   -11.9585 -11.9422 -11.9363 -11.9352 -11.9315 -11.9262 -11.9180 -11.9159
   -11.9103 -11.9067 -11.9062 -11.9057 -11.9048 -11.8994 -11.8959 -11.8932
   -11.8923 -11.8891 -11.8824 -11.8778 -11.8703 -11.8648  -7.9298  -7.6610
    -7.6124  -7.2183  -6.9747  -6.9668  -6.6880  -6.6054  -6.5897  -6.5510
    -6.4994  -6.3826  -6.3585  -6.3143  -6.2991  -6.2493  -6.2093  -6.1886
     0.0844   0.1032   0.3168   0.5798   0.8072   0.9918   3.8682   4.0410
     4.0751   4.2191   4.2410   4.2963   4.3693   4.4447   4.4568   4.4676
     4.5981   4.6159   4.6860   4.8044   4.8752   4.9984   5.0440   5.0939
     5.3674   5.6677   5.7680   5.7754   5.8993   5.9813   6.0052   6.3135
     6.3573   6.3897   6.4866   6.5433   6.6929   6.7892   6.8224   6.9036
     7.0890   7.3218   7.3511   7.3854   7.3935   7.4864   7.5094   7.6796
     7.7014   7.7789   7.9328   7.9926   8.0789   8.2021   8.2450   8.3066
     8.3295   8.4029   8.4062   8.4969   8.5329   8.6039   8.6821   8.8599
     8.9853   9.0892   9.3776   9.3870   9.4620   9.5350   9.6818   9.7596
     9.8797   9.8880   9.9521  11.7438  11.7621  11.7626  11.7975  11.8344
    11.8765  11.9409  11.9517  11.9780  12.4927  12.7702  12.8663  12.9062
    13.0424  13.1526  13.4715  13.9559  14.0039  14.0574  14.3466  14.4102
    14.4843  14.5149  14.9641  15.3032  15.3786  15.3906  15.4675  15.5442

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.0000 0.3849 0.1335 ( 11170 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8101 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0284 -12.0168 -12.0067 -11.9880 -11.9853 -11.9852 -11.9640 -11.9631
   -11.9486 -11.9475 -11.9353 -11.9339 -11.9291 -11.9273 -11.9206 -11.9130
   -11.9095 -11.9092 -11.9070 -11.9041 -11.9023 -11.9003 -11.8997 -11.8919
   -11.8885 -11.8871 -11.8866 -11.8757 -11.8708 -11.8631  -7.9003  -7.7790
    -7.4962  -7.1646  -7.1370  -6.8911  -6.7195  -6.6044  -6.5924  -6.4930
    -6.4521  -6.4269  -6.3487  -6.3377  -6.3331  -6.2110  -6.2057  -6.2054
     0.0510   0.2088   0.2188   0.6375   0.7635   1.0035   3.8920   3.9573
     4.1671   4.2020   4.2211   4.3466   4.3482   4.3756   4.4716   4.4820
     4.5782   4.6391   4.8003   4.8057   4.8444   4.9791   4.9853   5.1097
     5.3878   5.5805   5.6945   5.7758   5.8492   6.0737   6.1672   6.1851
     6.4118   6.4158   6.4764   6.5493   6.7215   6.7711   6.7968   6.9665
     7.1017   7.2380   7.2932   7.3725   7.4225   7.4942   7.5383   7.7065
     7.7698   7.8046   7.8711   7.8731   8.1367   8.2005   8.2126   8.3359
     8.3601   8.4123   8.4411   8.4624   8.4969   8.5684   8.6742   8.9355
     8.9557   9.1184   9.2690   9.4198   9.5290   9.5586   9.6538   9.7477
     9.8374   9.9390   9.9502  11.7353  11.7703  11.7872  11.7880  11.8114
    11.8855  11.9417  11.9558  11.9672  12.5314  12.7136  12.8169  12.9613
    13.0394  13.1573  13.5219  13.8739  13.9426  14.1196  14.1769  14.4113
    14.6624  14.7685  14.8014  15.2309  15.3138  15.4511  15.5016  15.5869

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333 0.5774 0.0000 ( 11157 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8101 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0015 -12.0014 -12.0014 -11.9977 -11.9977 -11.9976 -11.9520 -11.9520
   -11.9519 -11.9410 -11.9410 -11.9410 -11.9178 -11.9178 -11.9177 -11.9128
   -11.9128 -11.9127 -11.9034 -11.9033 -11.9033 -11.8999 -11.8999 -11.8999
   -11.8848 -11.8847 -11.8847 -11.8800 -11.8799 -11.8799  -7.4096  -7.4094
    -7.4088  -7.3150  -7.3143  -7.3142  -6.7070  -6.7069  -6.7068  -6.4854
    -6.4853  -6.4852  -6.3681  -6.3678  -6.3677  -6.2137  -6.2134  -6.2133
     0.3868   0.3868   0.3869   0.4415   0.4416   0.4416   3.9059   3.9059
     3.9064   4.3075   4.3076   4.3078   4.4772   4.4772   4.4775   4.8350
     4.8351   4.8352   4.9243   4.9244   4.9247   5.0553   5.0559   5.0559
     5.6952   5.6956   5.6956   5.9337   5.9340   5.9341   6.2650   6.2652
     6.2654   6.4699   6.4699   6.4700   6.8226   6.8226   6.8226   7.0551
     7.0552   7.0553   7.3279   7.3279   7.3281   7.6286   7.6289   7.6290
     7.7873   7.7874   7.7877   7.9717   7.9717   7.9720   8.3569   8.3569
     8.3572   8.5086   8.5087   8.5087   8.6103   8.6104   8.6106   9.1418
     9.1421   9.1422   9.2827   9.2827   9.2830   9.4649   9.4652   9.4653
     9.7929   9.7932   9.7934  11.8075  11.8077  11.8079  11.9080  11.9080
    11.9082  11.9848  11.9851  11.9853  12.6074  12.6078  12.6079  12.9889
    12.9890  12.9894  13.6159  13.6160  13.6160  13.8819  13.8819  13.8820
    14.7850  14.7851  14.7855  14.9893  14.9896  14.9900  15.4347  15.4358

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333 0.5774 0.1335 ( 11151 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8101 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0005 -12.0004 -12.0004 -11.9985 -11.9984 -11.9983 -11.9495 -11.9494
   -11.9493 -11.9440 -11.9440 -11.9439 -11.9156 -11.9156 -11.9155 -11.9136
   -11.9135 -11.9134 -11.9028 -11.9028 -11.9027 -11.9011 -11.9011 -11.9010
   -11.8832 -11.8831 -11.8831 -11.8808 -11.8808 -11.8808  -7.4036  -7.4035
    -7.4027  -7.3240  -7.3232  -7.3231  -6.6767  -6.6766  -6.6764  -6.5449
    -6.5448  -6.5447  -6.3202  -6.3199  -6.3197  -6.2297  -6.2292  -6.2291
     0.3998   0.3998   0.4000   0.4280   0.4282   0.4282   3.8963   3.8963
     3.8967   4.3592   4.3593   4.3595   4.4488   4.4489   4.4492   4.8167
     4.8168   4.8168   4.9351   4.9352   4.9354   5.0564   5.0570   5.0570
     5.6391   5.6393   5.6394   5.9624   5.9625   5.9627   6.2780   6.2780
     6.2782   6.5018   6.5019   6.5020   6.8349   6.8350   6.8352   7.0638
     7.0642   7.0642   7.2468   7.2468   7.2470   7.6329   7.6330   7.6332
     7.7584   7.7585   7.7585   8.0703   8.0704   8.0705   8.3701   8.3702
     8.3706   8.4812   8.4812   8.4813   8.5892   8.5893   8.5894   9.1543
     9.1545   9.1546   9.2468   9.2469   9.2471   9.4723   9.4725   9.4725
     9.8064   9.8067   9.8068  11.8088  11.8089  11.8092  11.8984  11.8985
    11.8985  11.9921  11.9924  11.9926  12.6106  12.6110  12.6111  12.9950
    12.9952  12.9956  13.6733  13.6734  13.6737  13.7680  13.7683  13.7684
    14.8818  14.8820  14.8825  14.9807  14.9807  14.9810  15.5000  15.5014

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k =-0.3333 0.1925 0.1335 ( 11170 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8100 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0311 -12.0047 -12.0033 -12.0006 -11.9968 -11.9720 -11.9692 -11.9542
   -11.9496 -11.9489 -11.9382 -11.9366 -11.9259 -11.9242 -11.9218 -11.9140
   -11.9105 -11.9097 -11.9061 -11.9060 -11.9041 -11.9005 -11.8961 -11.8929
   -11.8909 -11.8835 -11.8805 -11.8791 -11.8773 -11.8606  -7.8885  -7.8056
    -7.4246  -7.3066  -7.0546  -6.8746  -6.6973  -6.6473  -6.5552  -6.5254
    -6.4963  -6.3919  -6.3790  -6.3127  -6.2957  -6.2319  -6.2222  -6.1889
     0.0502   0.1356   0.3853   0.4657   0.9078   0.9385   3.8950   3.9666
     4.1335   4.1998   4.2882   4.3010   4.3453   4.4153   4.4380   4.5238
     4.5323   4.6653   4.7102   4.8048   4.8797   4.9790   4.9796   5.1256
     5.4707   5.5576   5.6846   5.7503   5.8475   6.0393   6.1776   6.2977
     6.2979   6.3895   6.4800   6.5505   6.7842   6.7989   6.8188   6.8694
     7.1536   7.1815   7.3620   7.3939   7.4587   7.4889   7.5752   7.5802
     7.7557   7.7740   7.7929   8.0342   8.0977   8.1780   8.2315   8.3338
     8.3468   8.3845   8.4784   8.4910   8.5364   8.6039   8.6580   8.8378
     8.9191   9.1320   9.2764   9.4689   9.5066   9.5742   9.6898   9.6954
     9.8810   9.9026   9.9519  11.7459  11.7506  11.7707  11.7907  11.8504
    11.8772  11.9331  11.9572  11.9709  12.5640  12.6541  12.8995  12.9414
    12.9702  13.2071  13.6239  13.7485  13.9707  14.1050  14.2692  14.3243
    14.5933  14.6421  15.1570  15.2359  15.2382  15.2788  15.3284  15.6812

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

          k = 0.3333-0.5774-0.1335 ( 11151 PWs)   bands (ev):

   -79.4735 -79.4735 -79.4731 -74.6885 -74.6880 -74.6877 -45.8702 -45.8702
   -45.8699 -45.8101 -45.8101 -45.8101 -45.8100 -45.8097 -45.8097 -41.2479
   -41.2474 -41.2472 -41.1986 -41.1986 -41.1982 -41.1982 -41.1979 -41.1979
   -12.0006 -12.0006 -12.0006 -11.9983 -11.9982 -11.9981 -11.9497 -11.9496
   -11.9495 -11.9437 -11.9437 -11.9437 -11.9159 -11.9159 -11.9158 -11.9133
   -11.9133 -11.9132 -11.9030 -11.9029 -11.9029 -11.9010 -11.9010 -11.9009
   -11.8832 -11.8831 -11.8830 -11.8809 -11.8809 -11.8808  -7.4029  -7.4028
    -7.4020  -7.3249  -7.3241  -7.3240  -6.6718  -6.6717  -6.6715  -6.5524
    -6.5523  -6.5521  -6.3147  -6.3143  -6.3141  -6.2324  -6.2319  -6.2317
     0.4002   0.4002   0.4004   0.4276   0.4277   0.4277   3.9005   3.9005
     3.9010   4.3561   4.3562   4.3564   4.4448   4.4449   4.4452   4.8200
     4.8200   4.8201   4.9291   4.9293   4.9294   5.0524   5.0529   5.0530
     5.6807   5.6809   5.6810   5.9290   5.9293   5.9294   6.2630   6.2632
     6.2634   6.5158   6.5158   6.5159   6.8218   6.8218   6.8218   7.1133
     7.1133   7.1136   7.2235   7.2237   7.2238   7.6051   7.6053   7.6055
     7.8240   7.8241   7.8241   8.0025   8.0026   8.0027   8.3359   8.3360
     8.3362   8.4826   8.4827   8.4827   8.6485   8.6486   8.6487   9.1593
     9.1596   9.1597   9.2572   9.2573   9.2576   9.4514   9.4517   9.4517
     9.8021   9.8024   9.8025  11.8084  11.8085  11.8087  11.9144  11.9145
    11.9147  11.9810  11.9814  11.9816  12.6331  12.6335  12.6336  12.9672
    12.9673  12.9677  13.6223  13.6224  13.6225  13.8225  13.8227  13.8228
    14.8793  14.8795  14.8796  14.9908  14.9910  14.9918  15.4818  15.4834

     occupation numbers 
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000
     1.0000   1.0000   1.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000
     0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000   0.0000

     the Fermi energy is    11.0897 ev

!    total energy              =   -3013.09242065 Ry
     Harris-Foulkes estimate   =   -3013.09242073 Ry
     estimated scf accuracy    <       0.00000028 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =   -1374.87329006 Ry
     hartree contribution      =     856.48884751 Ry
     xc contribution           =    -494.27606785 Ry
     ewald contribution        =   -2000.90100253 Ry
     Hubbard energy            =       0.46909227 Ry
     smearing contrib. (-TS)   =      -0.00000000 Ry

     total magnetization       =     0.00 Bohr mag/cell
     absolute magnetization    =    25.71 Bohr mag/cell

     convergence has been achieved in  15 iterations

     Writing output data file BiFeO3.save/

     init_run     :    113.27s CPU     24.53s WALL (       1 calls)
     electrons    :  17482.77s CPU   2938.43s WALL (       1 calls)

     Called by init_run:
     wfcinit      :     98.42s CPU     17.79s WALL (       1 calls)
     wfcinit:atom :      1.11s CPU      0.24s WALL (      16 calls)
     wfcinit:wfcr :     86.14s CPU     14.53s WALL (      16 calls)
     potinit      :      2.00s CPU      1.77s WALL (       1 calls)
     hinit0       :      4.60s CPU      3.04s WALL (       1 calls)

     Called by electrons:
     c_bands      :  16706.03s CPU   2801.47s WALL (      15 calls)
     sum_band     :    661.77s CPU    113.34s WALL (      15 calls)
     v_of_rho     :      4.44s CPU      0.79s WALL (      16 calls)
     v_h          :      1.32s CPU      0.24s WALL (      16 calls)
     v_xc         :      3.11s CPU      0.55s WALL (      16 calls)
     newd         :    115.37s CPU     22.67s WALL (      16 calls)
     mix_rho      :      2.46s CPU      1.09s WALL (      15 calls)

     Called by c_bands:
     init_us_2    :     19.09s CPU      3.27s WALL (     512 calls)
     ccgdiagg     :  15652.50s CPU   2624.11s WALL (     240 calls)
     wfcrot       :   1120.23s CPU    188.60s WALL (     240 calls)

     Called by sum_band:
     sum_band:bec :      0.76s CPU      0.12s WALL (     240 calls)
     addusdens    :    158.24s CPU     28.71s WALL (      15 calls)

     Called by *cgdiagg:
     h_psi        :  10706.68s CPU   1793.65s WALL (  100781 calls)
     s_psi        :   3190.44s CPU    534.99s WALL (  201338 calls)
     cdiaghg      :     23.85s CPU      4.77s WALL (     240 calls)

     Called by h_psi:
     h_psi:pot    :  10259.84s CPU   1717.85s WALL (  100781 calls)
     h_psi:calbec :   1666.76s CPU    279.48s WALL (  100781 calls)
     vloc_psi     :   3078.14s CPU    515.04s WALL (  100781 calls)
     add_vuspsi   :   5512.53s CPU    923.04s WALL (  100781 calls)
     vhpsi        :    430.97s CPU     72.98s WALL (  100781 calls)
     h_1psi       :  11316.30s CPU   1896.01s WALL (  100541 calls)

     General routines
     calbec       :   3349.80s CPU    562.84s WALL (  302599 calls)
     fft          :      4.92s CPU      1.86s WALL (     172 calls)
     ffts         :      0.57s CPU      0.17s WALL (      62 calls)
     fftw         :   3304.31s CPU    551.72s WALL (  328698 calls)
     interpolate  :      0.50s CPU      0.50s WALL (      32 calls)
     davcio       :      0.00s CPU      0.03s WALL (      16 calls)


     Hubbard U routines
     new_ns       :      8.88s CPU      1.49s WALL (      15 calls)
     vhpsi        :    430.97s CPU     72.98s WALL (  100781 calls)

     PWSCF        :     4h53m CPU    49m23.95s WALL


   This run was terminated on:  10:53:27  13Oct2018            

=------------------------------------------------------------------------------=
   JOB DONE.
=------------------------------------------------------------------------------=
